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PREFACE 


This publication is a collection of the presentations given at the NASA Computational 
Fluid Dynamics (CFD) Conference held at NASA Ames Research Center, Moffett Field, 
California, March 7-9, 1989. The objectives of the conference were to disseminate CFD 
research results to industry and university CFD researchers, to promote synergy among 
NASA CFD researchers, and to permit feedback from researchers outside NASA on 
issues pacing the discipline of CFD. The focus of the conference was on the application 
of CFD technology but also included fundamental activities. The conference was 
sponsored by the Aerodynamics Division, Office of Aeronautics and Space Technology 
(OAST), NASA Headquarters, Washington, DC 20546. 

The conference consisted of twelve sessions of papers representative of CFD research 
conducted within NASA and three non-NASA panel sessions. For each panel session, the 
panel membership consisted of industry and university CFD researchers. A summary of 
the comments made during the panel sessions have been included in this publication. 

The conference proceedings are published in two volumes. Volume 1 contains the papers 
presented in Sessions I- VI; Volume 2 contains those given in Sessions VII-XII. Each 
volume contains the same front matter, and each contains a list of attendees as an 
appendix. 
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NASA CFD Conference 
NASA Ames Research Center 
March 7-9, 1989 

Panel Sessions Summary 

The NASA CFD Conference was held at Ames Research Center on March 7-9, 
1989. To conclude each day's presentations, a panel session with participation 
from the audience furnished a great deal of excellent feedback from the industry 
and academic communities. During the conference it was evident that the panel 
members proffered comments only after having spent considerable time in 
preparing them. 

The members of the panel sessions are listed below: 

March 7 P. Rubbert - Boeing Commercial Airplanes 

R. Melnik - Grumman Aerospace Corporation 
D. Whitfield - Mississippi State University 

March 8 I. Bhateley - General Dynamics - Fort Worth Division 

R. Agarwal - McDonnell Douglas Research Laboratories 
R. MacCormack - Stanford University 

March 9 V. Shankar - Rockwell International Science Center 
J. Carter - United Technologies Research Center 
A. Jameson - Princeton University 

The crucial comments from the three panel sessions have been combined and are 
summarized as follows: 

- NASA's CFD program is now too heavily focused on applications: program 
balance has swung from fundamentals (1970's) to applications (1980's) 

- Three critical "needs" emerged: 

(1) More algorithm research is needed; especially for Navier-Stokes 
solvers with unstructured grids 

(2) More research is required on geometric modelling; need rapid, 
accurate, and effective surface definition techniques 

(3) More research is needed on grid generation methods with the focus 
on speed, efficiency, and grid quality to reduce set up time and 
complexity 

- Developers of CFD need to understand the needs of the users; designers of 
aerospace vehicles have requirements that are different than the CFD 
researchers perceptions 

- Industry needs more reliable and cost effective CFD tools 
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Additional detail comments from the three panel sessions are listed below: 

- CFD has matured during the last decade and is being used to solve real 
problems; however, industry lacks confidence in Navier-Stokes solutions 

- Industry needs codes that have been validated to increase confidence in 
CFD technology 

- Improved communality between codes would increase usability; standards 
are needed 

- Improved data storage, networking, data transfer, and graphics required to 
assimilate information provided by CFD 

- Improved turbulence modeling for separated flows 

- Accurate prediction of drag for complete powered aerospace vehicles 

- Develop multidisciplinary CFD technology with optimization capability 

- NASA must maintain focus on technology development and high risk 
research 

- Technology transfer is not complete until design engineers are using CFD 
codes successfully 

- Industry needs NASA to improve CFD technology for codes simpler than 
Navier-Stokes solvers 

- NAS program has been extremely helpful to industry in transferring CFD 
technology 

- Industry needs to be more aggressive in their use of CFD 

- Improved understanding of CFD by design engineers required; cooperative 
programs or workshops were suggested to bring CFD researchers and 
designers together 

- Design cycle time needs to be reduced with CFD; codes must be cost 
effective, reliable, and useable, and robust to work at flight Reynolds 
Numbers 

- Improved coordination/reduced overlap of CFD applications between NASA 
centers 
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A COMPARATIVE STUDY OF NAVIER-STOKES CODES 
FOR HIGH-SPEED FLOWS 


David H. Rudy, James L. Thomas, Ajay Kumar, Peter Gnoffo 
NASA Langley Research Center 

and 

Sukumar R. Chakravarthy 
Rockwell International Science Center 


A comparative study has been made with four different codes for solving the compressible 
Navier-Stokes equations using three different test problems. The first of these cases was hypersonic 
flow through the P8 inlet, which represents inlet configurations typical of a hypersonic airbreathing 
vehicle. The free-stream Mach number in this case was 7.4. This 2-D inlet was designed to provide 
an internal compression ratio of 8. Initial calculations were made using two state-of-the-art finite- 
volume upwind codes, CFL3D and USA-PG2, as well as NASCREN, a code which uses the unsplit 
finite-difference technique of MacCormack. All of these codes used the same algebraic eddy-viscosity 
turbulence model. In the experiment, the cowl lip was slightly blunted; however, for the computations, 
a sharp cowl leading edge was used to simplify the construction of the grid. Although the overall 
features of the inlet flow field were predicted reasonably well, discrepancies between the computed and 
experimental profiles of pitot pressure and total temperature were found for all of the codes. Calcula- 
tions were then made with the two finite-volume upwind codes using a patched-grid approach which 
allowed the use of a blunt cowl tip. This produced a more accurate location for the cowl leading-edge 
shock and thus an improvement in the location of the resultant shock reflections within the inlet. 

The second test problem was the supersonic (Mach 3.0) flow in a three-dimensional comer 
formed by the intersection of two wedges with equal wedge angles of 9.48 degrees. The flow in such 
a comer is representative of the flow in the comers of a scramjet inlet. Calculations were made for 
both laminar and turbulent flow and compared with experimental data. The three-dimensional versions 
of the three codes used for the inlet study (CFL3D, USA-PG3, and SCRAMIN, respectively) were 
used for this case. For the laminar comer flow, a fourth code, LAURA, which also uses recently- 
developed upwind technology, was also utilized. It was found that the complex flow structure was 
qualitatively predicted by all of the codes. Furthermore, with a sufficiently-refined grid, all four codes 
gave identical results for the laminar case. 

The final test case is the two-dimensional hypersonic flow over a compression ramp. In this 
case, the flow is laminar with a free-stream Mach number of 14.1. In the experiment, the ramp angle 
was varied to change the strength of the ramp shock and the extent of the viscous-inviscid interaction. 
Calculations have been made for the 24-degree ramp configuration which produces a large separated- 
flow region that extends upstream of the comer. All of the codes predicted the strong-interaction struc- 
ture observed experimentally. However, a wide variation in the predicted extent of separation was 
found for the grid that was used. 
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DESCRIPTION OF CODES 


CFL3D 

• Written by James L. Thomas 

• Upwind-biased Roe scheme for pressure 
and convective terms 

• Central-difference treatment of viscous terms 

• Approximate-factorization time differencing 


USA-PG3 

• Written by Sukumar Chakravarthy 

• Roe scheme for convective and pressure 
terms 

• Central-difference treatment of viscous terms 

• Approximate-factorization time differencing 


DESCRIPTION OF CODES 


LAURA 

® Written by Peter A. Gnoffo 

• Symmetric TVD scheme with Roe averaging 

• Alternating directional-sweep Gauss-Seidel 
relaxation 

• Laminar flow only at present time 

NASCRIN and SCRAMIN 

• Written by Ajay Kumar 

• Explicit unsplit MacCormack scheme 


TEST PROBLEMS 


P8 Inlet 

• 2-D, transitional 

• Experimental data: Gnos et al. 


Symmetric-wedge corner 

• 3-D, laminar and turbulent 

• Experimental data: West and Korkegi 


Hypersonic compression corner 

• 2-D, laminar 

• Experimental data: Holden and Moselle 


P8 INLET 


Wedge shock 


Cowl shock 



Throat 


GRIDS NEAR INLET ENTRANCE 


Patched-grid computations 



COMPUTED INLET PRESSURE 


USA-PG2 

Transitional computation 
Sharp cowl leading edge 



CENTERBODY PRESSURE 


Blunt cowl leading edge 



PITOT PRESSURE PROFILE 


Y w (in.) 


x = 40.9 in. 

Blunt cowl leading edge 



CORNER FLOW 


IVL=3.0 



1 Corner shock 

2 Internal shock 

3 Slip line 

4 Wall shock 



LAMINAR CORNER FLOW 


x = 0.0733 m. 



Density contours 
USA-PG3 


O Experiment 



Sidewall pressure 


LAMINAR CORNER FLOW 

Sidewall pressure 
x = 0.0733 m. 


Conical flow 


121 x 121 grid 



x 


HYPERSONIC COMPRESSION CORNER 


Streamlines 
24-degree ramp 

USA-PG2 


SURFACE PRESSURE 


24° ramp 
201 x 60 grid 
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SKIN FRICTION 


24° ramp 
201 x 60 grid 
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SUMMARY AND CONCLUSIONS 


A comparative study of four computer codes for 
the Navier-Stokes equations was performed 
using three test problems for high-speed flow. 


For the P8 inlet, good agreement was found 
between the codes. Differences between 
computation and experiment were possibly due 
to experimental effects not included. 

Flow structure in 3-D corner was predicted well 
by all of the codes. Good agreement was 
found between the solutions for all four codes. 


For the hypersonic compression ramp, the 
codes predicted the strong interaction structure 
observed experimentally. A wide variation in the 
predicted extent of separation was found. 
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MODELING OF HIGH SPEED CHEMICALLY REACTING FLOW-FIELDS 

J. P. Drummond, M. H. Carpenter* and H. Kamath* 

NASA Langley Research Center 
Hampton, VA 23665 

ABSTRACT 

The SPARK3D and SPARK3D-PNS computer programs have been developed to model 
3-D supersonic, chemically reacting flow-fields. The SPARK3D code is a full Navier-Stokes 
solver, and is suitable for use in scramjet combustors and other regions where recirculation 
may be present. The SPARK3D-PNS is a parabolized Navier-Stokes solver and provides 
an efficient means of calculating steady-state combustor far-fields and nozzles. Each code 
has a generalized chemistry package, making modeling of any chemically reacting flow 
possible. 

Research activities by the Langley group range from addressing fundamental theoretical 
issues to simulating problems of practical importance. Algorithmic development includes 
work on higher order and upwind spatial difference schemes. Direct numerical simulations 
employ these algorithms to address the fundamental issues of flow stability and transition, 
and the chemical reaction of supersonic mixing layers and jets. It is believed that this 
work will lend greater insight into phenomenological model development for simulating 
supersonic chemically reacting flows in practical combustors. Currently, the SPARK3D and 
SPARK3D-PNS codes are used to study problems of engineering interest, including various 
injector designs and 3-D combustor-nozzle configurations. Examples, which demonstrate 
the capabilities of each code are presented. 
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OVERVIEW 


GROUP OBJECTIVES 
DESCRIPTION OF 2 AND 3-D CODES 
3-D FULL NAVIER-STOKES 
3-D PARABOLIZED NAVIER-STOKES 
2 AND 3-D MIXING ENHANCEMENT 


CONCLUSIONS AND DIRECTIONS 


OBJECTIVES 


THEORETICAL ISSUES 


- NUMERICS FOR SUPERSONIC CHEMICALLY 
REACTING FLOWS 

- Numerical Efficiency 

- Numerical Accuracy 

- Robustness 


- PHYSICAL ISSUES 

- Mixing Enhancement and Combustion 

- Transition to Turbulence 

- Phenomenological Turbulence Models 


APPLICATIONS 

- 3-D COMBUSTORS AND NOZZLES 

- SPARK3D 

- SPARK3D-PNS 


SPARK 3-D 


3-D NAVffiR-STOKES AND CHEMISTRY 
FULLY ELLIPTIC STRUCTURE 
VISCOUS OR INVISCID CAPABILITIES 
TIME-ACCURATE OR LOCAL RELAXATION 
GENERALIZED CHEMISTRY ROUTINES 


LOW STORAGE FORMAT 


GENERALIZED CHEMISTRY ROUTINE 


REAL GAS THERMODYNAMICS MODELS 

- POLYNOMIAL FITS FOR C p , C v , S, G 

- KINETIC THEORY BASED DIFFUSION 
MODELS 

- Sutherlands Law for jij and kj 

- Wilkes law for |i and k 

- Binary or multicomponent diffusion models 


CHEMISTRY MODELS 

- FROZEN 

- EQUILIBRIUM 

- FINITE RATE CHEMISTRY (NS=9, NR=18) 


NUMERICAL METHODS 


TRANSFORMED COORDINATES (0T|,0 

- GEOMETRIC CONSERVATION LAW TERMS 
INCLUDED 


TEMPORAL INTEGRATION (2ND ORDER) 

- EXPLICIT FORMULATION FOR HYDRO- 
DYNAMIC TERMS 

- Allows local time stepping 

- EXPLICIT/IMPLICIT FORMULATION FOR 
CHEMICAL SOURCE TERMS 


FINITE DIFFERENCE SPATIAL DISCRETIZATION 

- STANDARD MACCORMACK (2ND ORDER) 

- GOTTLIEB MACCORMACK (4TH ORDER) 

- COMPACT MACCORMACK (4TH ORDER AT 
S.S.) 

- UPWIND (3RD ORDER) 


CORNER FLOW 



1 Corner shock 

2 Internal shock 

3 Slip line 

4 Wall shock 



Figure 3.- Symmetric- wedge corner and schematic 

of corner flow. 


CORNER FLOW 

Symmetric- Wedge Corner 



Density Contours (Laminar Flow) 


REARWARD STEP 


z 



L = 7.0 cm, H = 1 .8 cm, W = 2.9 cm 
L| = 0.7 cm, l_2= 1 .2 cm, S = 0.3 cm 


SPARK3D 

REARWARD STEP COMPARISON 

Mach = 2.05 
101 x 41 x 25 GRID 

Qr = 0.39 


JET PENETRATION IN X-Z PLANE AT Y = 



JET SPREAD IN X-Y PLANE AT Z = D 
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EXPERIMENTAL COMPARISON 

Jet Penetration 
61 x 41 x 25 GRID 
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EXPERIMENTAL COMPARISON 

Jet Penetration 


61 x 41 x 25 GRID 
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REARWARD STEP COMPARISON 

Jet Spread 
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SPARK3D-PNS 


EXTENSION OF SPARK3D 


EFFICIENT SOLUTIONS OF STEADY 3-D PNS 
EQUATIONS 


FOR USE IN THE COMBUSTOR FAR-FIELD AND 
NOZZLE 


INTEGRATION SCHEME BASED ON 2ND ORDER 
MACCORMACK ALGORITHM 


u> 

u> 


CORNER FLOW 



BURROWS KURKOV 


It involves Mach 1 H 2 injection in a Mach 2.44 vitiated-air stream. The flow 
conditions correspond to the experiments of Burrows and Kurkov. For the 
H 2 jet: u = 1216 m/s, T = 254 K and /# 2 = 1.0. For the vitiated-air stream: 
u = 1764 m/s, T = 1270 K, fo 2 = 0.258, /n 2 = 0.486 and fn 2 o = 0.256. 


Test section initial 
measurement station, 
0 cm 


Test section intermediate 
measurement station, 

18i3cm Test section exit 

measurement station, 
35. 6 cm 
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injection 

Static pressure ports 
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MIXING ENHANCEMENT STUDIES 


FUEL-AIR MIXING DECREASES WITH INCREAS- 
ING MACH NUMBER 


MIXING ENHANCEMENT MECHANISMS ARE 
REQUIRED AT HIGH COMBUSTOR MACH 
NUMBER 

- PLANAR SHOCKS \ EXCITATION 

- CURVED SHOCKS I OF 

- SWIRL [ UNSTABLE 

- ACOUSTIC FORCING J MODES 


ENHANCEMENT BY SHOCKS IS EXAMINED IN 
THIS STUDY 


ORIGINAL PAGE IS 

OF POOR QUALITY 


Cj 





STRUT MODIFICATION FOR IMPROVED COMBUSTION EFFICIENCY 


Conventional strut 


Modified strut 


Water mass fractio 


Conventional strut 


Water mass fraction 


Modified strut 


NASA 


L-88-9164 




COMBUSTION EFFICIENCY 


STRUT WITH TRANSVERSE - PARALLEL JETS 



MODEL COMBUSTOR 


i s i 



w 


L = 20.0 cm, W = 1 0.0 cm 
S = 2.0 cm, L 1 = 1 .5 cm, D = 3.5 mm 

T = 1 000 K, P = 0.5 atm, V = 1 500 m/s 


CONCLUSIONS 


COMPUTER PROGRAMS DEVELOPED TO 
MODEL 3-D SUPERSONIC CHEMICALLY 
REACTING FLOWFIELDS. 


FAVORABLE VALIDATION AGAINST AVAIL- 
ABLE EXPERIMENTAL RESULTS 


BEING USED EXTENSIVELY IN THEORETICAL 
STUDIES AND IN ENGINEERING ENVIRON- 
MENTS 


DIRECTIONS 


ALGORITHMS 

- HIGHER ORDER AND COMPACT 

- UPWINDING 

- PARALLEL COMPUTING 


TRANSITION - TURBULENCE MODELING 


- COMPARISON WITH LINEAR STABILITY 
CODES 

- DIRECT SIMULATIONS OF SUPERSONIC 
JETS AND MIXING LAYERS 

- "DATABASE" FOR PHEN OMIN OLOGIC AL 
MODELS 


PRODUCTION CODE SUPPORT 
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Three-Dimensional Calculation of Supersonic 
Reacting Flows Using an LU Scheme 

Sheng-Tao Yu 
Y-L Peter Tsai 
Jian-Shun Shuen 

Sverdrup Technology, Inc . 

NASA Lewis Research Center 
Cleveland, Ohio 


ABSTRACT 

A new three-dimensional numerical program incorporated with 
comprehensive real gas property models has been developed to sim- 
ulate supersonic reacting flows. The code employs an implicit finite 
volume, Lower-Upper (LU) time-marching method to solve the com- 
plete Navier-Stokes and species equations in a fully-coupled and very 
efficient manner. A chemistry model with nine species and eighteen 
reaction steps are adopted in the program to represent the chemical 
reaction of H2 and air. To demonstrate the capability of the program, 
flow fields of underexpanded hydrogen jets transversely injected into 
supersonic air stream inside the combustors of scramjets are calcu- 
lated. Results clearly depict the flow characteristics, including the 
shock structure, separated flow regions around the injector, and the 
distribution of the combustion products. 
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PREVIOUS WORK 


Develop fluid dynamic code to study the mixing and chemical 
reactions inside the combustors of ramjets. 

A two-dimensional computer code (RPLUS) using LU scheme 
to simulate chemical reacting flows has been developed: 

— Fully coupled and very efficient. 

— Finite rate chemistry. 

— A fast equilibrium chemistry package. 

Previous work done by using RPLUS 2D code: 

— Hypersonic inlet flows at Mach 10 and 13. 

— The combustion of a hydrogen jet transversely injected into 
a supersonic air stream. 

— Hydrogen air mixture passing a 10° ramp. 


ANALYSIS 


• 3D Navier-Stokes and N a -1 species equations. 


• Brabbs’ finite rate chemistry model. 

9 species: H 2 , H, OH, H 2 0, O, 0 2 , H0 2 , H 2 0 2 , and N 2 . 


- 18 reaction steps. 

• Implicit treatment of source terms in species equations. 

• No source term in energy equation for either exothermic or 
endothermic reactions. 

• Temperature and pressure are calculated iteratively from the 
following equations: 



N. 


t'=l 



T 
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THERMODYNAMIC AND TRANSPORT PROPERTIES 


• C p , thermal conductivity, and viscosity of individual species 
are obtained from fourth order polynomial of T. 

• Binary mass diffusivity is calculated using the Chapman- 
Enskog theory with Lennard-Jones intermolecular potential en- 
ergy functions. 

ON 

• Wilke’s mixing rule is used for the transport properties of gas 
mixture. 

• Mass concentration weighing is used for the C p of gas mixture. 

TURBULENCE MODEL 

• Baldwin-Lomax model. 


Sc t = Prt = 0.9 


NUMERICAL SCHEME 


• Time-marching LU scheme. 

• Implicit treatment of source terms. 
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NUMERICAL SCHEME 

• Two operators for 3D calculation. 

• Operation count comparable to explicit schemes. 


NAQijk - AtTijk AQijk + 
At 


Ax 

At 
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PROGRAM VECTORIZATION 


• LHS can be vectorized on planes normal to the sweeping direc- 
tion. 



Lower-Sweep 


Upper-Sweep 
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CONVERGENCE HISTORY 


• Coarse grid (30x25x30) 



• Present calculation: 

- Residuals reduced by 4 orders of magnitude. 

- Single injection case: (61x39x43) 8 MW, 5 hrs on Cray 2 
Dual Injection case: (81x39x43) 11 MW, 6 hrs on Cray 2 
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GEOMETRY AND FLOW CONDITION 
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MASS-FRACTION CONTOUR [Q = 1.02,60 x 100 x 50) 


Jet 



J. C. MacDaniel 


MASS-FRACTION CONTOUR ( Two Ports. Q = 1.02, 70 x 100 x 50) 
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FLOW CONFIGURATION 



Mach Number 



x-0 .615cm 


x=l . 16cm 


x=1.61cm 


x=2. 19cm 
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CONCLUDING REMARKS 


A new three-dimensional computer code for high speed reacting 
flows has been successfully developed. 

The code is very efficient because the LHS can be vectorized 
on the planes normal to the sweeping directions. 

The validity of the code is assessed by the comparison between 
the calculated results and MacDaniel’s experimental data for 
non-reacting mixing flows. 

The code is applied to chemically reacting flows of H2 jets trans- 
versely injected into hot airstream. 

- Results clearly depict the shock structure, recirculations, 
and species distribution due to chemical reaction. 

- Two-hole injection allows deeper penetration of the fuel jet 
and more complete combustion of H2 and air. 
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Progress in Computing 
Nozzle/Plume Flow Fields 

Stephen M. Ruffin 

NASA Ames Research Center, Moffett Field, CA 

Ethiraj Venkatapathy 

Eloret Institute, Sunnyvale, CA 

William J. Feiereisen 

NASA Ames Research Center, Moffett Field, CA 

Seung-Ho Lee 

Eloret Institute, Sunnyvale, CA 


Abstract 

The long-term goal of this work is to develop the capability to predict 
chemically-reacting, multi-stream nozzle and plume flow fields. Two basic 
Navier-Stokes solvers, including the widely used F-3D code, are upgraded 
to include several upwind difference schemes and portable chemistry packages. 
Current computational capabilities for solving eqilibrium single-stream and 
multi-stream, frozen gas and finite rate chemistry problems are described. 

A variety of complex nozzle and plume flows have been computed. Solutions 
presented herein include axisymmetric plume flow for ideal and equilibrium air, 
3-D NASP nozzle /afterbody flow, and an internal nozzle calculation comparing 
various finite-rate chemistry packages. 
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Motivations, Objective, and Applications 


Motivations 

• Flow in nozzle and propulsive plume strongly 
influence performance of hypersonic vehicles. 

• Ground-based facilities can not fully simulate these 
flow fields. 

• Validated, accurate and efficient flow codes are needed 
to design integrated propulsion systems of hypersonic 
vehicles. 


Objective 

• Develop computational ability to predict chemically 
reacting, multi-stream, nozzle/plume flow fields. 

Applications 

• NASP, Space Shuttle, general plume flows and 
signatures 

/ 


SCHEMATIC OF AXISYMMETRIC PLUME FLOWS 



Plume Flow in Subsonic External Stream 



Plume Flow in Supersonic External Stream 


— 
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Nozzle/Plume Flow Solvers 

Difference 

Scheme 

2-D Axisym. Eqns. 
Flux-Split in 2-D 

3-D Eqns. 
Flux-Split in 1-D 
Central Diff. in 2-D 

3-D Eqns. 

Roe’s upwind 
averaging in 3-D 

Implicit Soln. 
Algorithm 

Beam & Warming in 
"delta form" 

Two-Factor Flux-Split 

LDU-ADI 

Viscous Terms 

in 1-D 

in 1-D 

in 3-D 

Chemistry 

Ideal Gas 

Equil. single-stream 
Equil. multi-stream 
Frozen 

Ideal Gas 
Finite-Rate 

Ideal Gas 

Equil. single-stream 

Calibration 

Ideal Gas: Mach 
disk location for 
jet plumes 
Equil. 1-stream: 

Shock standoff dist., 
and shape and species 
distribution for blunt 
body flow. 

Ideal Gas: Numerous 
flows including Space 
Shuttle, Mach No., 
pressure and Mach 
disk location for plumes 
Finite-Rate: H 2 -air 
species distribution for 
axsymm. nozzle. 

Ideal Gas: Numerous 
flows including Space 
Shuttle, Mach No., 
pressure and Mach 
disk location for plumes 
plume structure in 
quiescent freestream 




J 


O' 

u> 


Nozzle/Plume Flow Solvers (cont.) 


• All calculations performed for laminar flow 

• Adaptive grid routine used for all 2-D cases 

(3-D adaptive grid routine recently completed) 


• 3-D code can handle multiple, patched-grid regions 
and generalized specification of boundary conditions 
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Equilibrium and Frozen Gas Chemistry 


Constant Elemental Composition (i.e.. single-stream) 
• Uniform gas throughout flow field 

Equilibrium Gas 

- table look-up, curve fits 

- relatively fast 

- must be pre-calculated 


Variable Elemental Composition (i.e.. multi-stream) 

• Solution of species partial density equations allows for 
convection, diffusion and mixing of multiple streams 

Equilibrium Gas 

- Gordon-McBride free energy minimization 

- arbitrary gas mixtures 
Frozen Chemistry 

- No reactions 
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Finite Rate Chemistry 


• H 2 - air chemistry currently being used 

- JANNAF plume model (Dash and Pergament, 1980) 

- 7 Species: H, H 2 , H 2 O, O, OH, 0 2 , N 2 

- 8 Reactions 

• Specific heats from curve fits (Gordon and McBride data) 

• Binary diffusion 

- Fick's Law with D s j given by Reid, et.al., 1977) 

• Non-catalytic wall 

• Chemistry Coupling 

- Fully-Coupled 

Invert large matrix of order 5 + (N SDGCi6S - 1) per node 

- Loosely-Coupled _ P 

Locally constant Y 

P not affected by species concentration variations 

Invert two smaller matrices of order 5 and N spe cies - 1 per node. 
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Types of Nozzle/Plume Flows Computed 



Main complexity is due to: 



Fluid Dyn. 

Body 

Flow 


Cases 

Features 

Geom. 

Chem. 

References 

1. Axisymmetric Single-Nozzle Plume 





A. Not including afterbody 





Ideal gas, sub. freestream 

X 



AIAA 89-0129 

Ideal gas, sup. freestream 

X 



AIAA 88-3158 

B. Including afterbody 





Ideal Gas 

X 

X 


AIAA 88-2636 

Equil. single-stream 

X 

X 

X 

AIAA 88-2636 

Equil. multi-stream 

X 

X 

X 

in preparation 

Frozen 

X 

X 

X 

in preparation 

II. 3-D Multiple-Nozzle Plume 





A. Not including afterbody 

X 



AIAA 88-3158 

B. Including afterbody 

X 

X 


AIAA 89-0129 

III. Hypersonic Nozzle/Afterbody 





A. 2-D Geometry 





Subsonic freestream 

X 



AIAA 89-0446 

Supersonic freestream 

X 



AIAA 89-0446 

B. 3-D Geometry 





Plume region alone 

X 

X 


AIAA 89-0446 

Complete geometry ◄- 

X 

X 


AIAA 89-0446 

JV. Internal Nozzle, finite rate chem. 



X 

in preparation j 


Axisymmetric Plume Flow in Quiscent Air 
Pj/Pinf = 21.9 - Mj = 1.5 


Upper Half - Computed Mach Contours with Adapted Grid 



Lower Half - Experimetnal Shadowgraph - (NASA TR R 6) 


Axlsymmetric Plume Solutions 
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A Flume and Body 
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Comparison of Shock Reflection Location vs Altitude 




PRECEDING PAGE BLANK NOT FILMED ORi'QiNAL PAGE IS 

OF POOR QUALITY 




PAGE BLANK NOT FILMED 


NOZZM'i GEOMETRY 
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Concluding Remarks 


• A variety of nozzle/plume flows have been successfully 
computed. 

- Complex flow structures and geometries have been studied 

- Good agreement with experimental and other computed 
results has been found for plume structure and Mach disk 
location. 

• Equilibrium single-stream, equil. multi-stream, finite rate, 
and frozen chemistry packages have been incorporated. 

• Turbulence modeling for complex plume flows is currently 
being studied 

• 3-D code will be validated against hypersonic nozzle/afterbody 
experimental data 

• Significant progress has been made toward goal of a single code 
capable of solving complex chemically reacting nozzle/plume 
flow fields. 
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APPLICATION OF CFD CODES FOR THE SIMULATION OF 
SCRAMJET COMBUSTOR FLOWFIELDS 


Tawit Chitsomboom 

NASA Langley Research Center/Vigyan Research Associates, Inc. 
Hampton, Virginia 23665-5225 

and 


G. Burton Northam 
NASA Langley Research Center 
Hampton, Virginia 23665-5225 


Office of Aeronautics and Space Technology (OAST) 

NASA Computational Fluid Dynamics Conference 
Mountain View, California 
March 7-9, 1989 

ABSTRACT 

An overview of CFD activities in the Hypersonic Propulsion Branch is 
given. Elliptic and PNS codes that are being used for the simulation of 
hydrogen-air combusting flowfields for scramjet applications are discussed. 
Results of the computer codes are shown in comparison with those of the 
experiments where applicable. Two classes of experiments will be presented: 
(a) parallel injection of hydrogen into vitiated supersonic air flow; and (b) 
normal injection of hydrogen into supersonic crossflow of vitiated air. 
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INTRODUCTION 

A number of CFD codes are being developed and 
validated at the Hypersonic Propulsion Branch, 

Fluid Mechanics Division, NASA Langley Research Center 

These CFD codes are being applied to solving flowfields 
inside hydrogen-fueled scramjet combustors 

Fully elliptic codes 

• 2D, Axisymmetric, 3D 

• Global two-step finite rate combustion model 

« Use unsplit explicit MacCormack Algorithm with point-implicit 
chemistry source terms 

PNS Codes 

• 2D, 3D 

• Global two-step finite rate combustion model 

• Implicit, coupled space marching procedure 


3D-ELLIPTIC CFD CODE 


• Unsplit explicit MacCormack finite-difference algorithm 

• Solve 9 PDE's by finite-difference technique 

• 1 Density 

• 3 Momentum 

• 1 Energy 

• 4 Species (O 2 , H 2 , H 2 O, OH) 


• Treat chemical source terms implicitly to alleviate 
stiffness associated with fast chemistry 

• 2D-Layer design for efficient use of memory 


COMBUSTION AND TURBULENCE MODELS 


# Four-species two-step finite-rate model (Rogers-Chinitz) 

kfl 

Ho + Oo ^ 20H 
* kbl 

kf2 

20H + Hp ^ 2HpO 

kb2 


kf, kb = kf, kb (<|> ,T) 

• Chemical source terms are obtained by applying the 
law of mass action 

• Prandtl mixing length hypothesis for the jet mixing process 

jit = kpcoL 2 
k » 0.02 


• Baldwin-Lomax model for near-wall turbulence 


RESULTS OF 2D ELLIPTIC CODE 


• Supersonic coaxial burner of hydrogen-air 

• Hot vitiated outer air stream co-flowing with cold 
inner hydrogen fuel, exiting into a quiescent air 

• Use "CARS" technique to obtain data point for 
oxygen, nitrogen and temperature 

• Compare CARS data with CFD results 


Water out 





CENTERLINE DISTRIBUTION OF O2 AND N 2 SPECIES 

Symbol : CARS data 
Lines : CFD data 
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0 2 Density 




X-inch 

• Wavy distributions due to shock waves reflection 
from the free stream 


• Flame was, infact, ignited by shock wave 
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RESULTS OF 2D ELLIPTIC CODE (CON’T) 

Comparison of radial distributions at 
various streamwise stations 

o CARS CFD 




(Temperature, K) 


3D-C0DE VALIDATION EXPERIMENT 

• Subscale scramjet combustor (3.5"W, 4"L, 1 .8"H) with : 

• Backward facing steps 

• Normal injections of H 2 pilot fuel, Silane ignitor fuel 
and H 2 primary fuel 

• 2° , Nozzle (4 ft long) 

• Inflow gas was Mach 2 vitiated air at high enthalpy. 

• Static pressures were measured at various points 
along the surfaces 
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CONTOURS AND VECTORS PLOT THROUGH 
THE MAIN JET: NO COMBUSTION 

Flow conditions: T s = 81 OK, P s = 1 atm, M = 2;<(> = 0.43 



Pressure 

contours 




Velocity 

vectors 




Hydrogen 

contours 


PRESSURE DISTRIBUTIONS ALONG BOTTOM WALL 

AT MIDPLANE 

Flow conditions: T s = 1200K, P s = 1 atm, M = 2; <t> = 0.48 

O Experiment, reacting 
O Experiment, non-reacting, no jets 



X/Gap 


• Mixing and combustion efficiency at exit was 
underpredicted by 20% 

• Prandtl mixing length hypothesis seems to work well 


HYDROGEN CONTOURS AND VELOCITY VECTORS 
AT X = 1.5" DOWNSTREAM OF MAIN JET 

Flow conditions: T s = 1200K, P s = 1 atm, M = 2; <(> = 0.48 




Kidney shape profile 
Counter rotating vortices 


3D-PNS CODE 


Applicable for supersonic high Reynolds number flow 

More efficient than full Navier-Stokes procedure 

Space marching technique using steady version of 
beam-warming algorithm (Vigneron, et al; Schiff-Steger) 

4 Species, 2-step combustion model (Rogers-Chinitz) 
for hydrogen-air 

Fluid and chemical species are solved together 
creating block (9 x 9) tridiagonal system 

Baldwin-Lomax turbulence model 

See also AIAA Paper 84-0438 


SCHEMATIC OF THE MODEL PROBLEM 
FOR THE 3D-PNS CODE 



Inflow conditions: Mach = 4, T = 900K, 
P = 1 atm, Equivalence ratio = 1 


FLOW FIELD AT THE MIDPLANE 
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Velocity vectors 



CONCLUDING REMARKS 


Computer codes are being developed and validated 
at HPB, FMD, NASA-LaRC 

Results so far have shown that MacCormack algorithm, 
despite some drawbacks, can give good 
physical results 

Simple Prandtl mixing length scheme is both viable 
and economical 

More sophisticated combustion and turbulence 
models can be used at the expense of computer time 
which is very costly for any 3D-reacting flow simulation 

3D-PNS code is an economical alternative to the 
costly elliptic code provided that certain 
conditions are met 

Need to make the PNS algorithm more robust 
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HYPERSONIC CFD APPLICATIONS AT NASA LANGLEY 
USING CFL3D AND CFL3DE 


Pamela F. Richardson 
NASA Langley Research Center 
Hampton, Virginia 
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CFL3D/CFL3DE 


• Time-dependent conservation law form of compressible Euler and 
Navier-Stokes equations 

Upwind-biased spatial differencing (Flux Vector Splitting - FVS and Flux 
Difference Splitting - FDS) 

• Thin-layer, finite-volume implementation with algebraic turbulence 
model 

Zonal grids - longitudinally patched (for hypersonic flows) 

- CFL3D 

• 3-factor implicit time advancement algorithm 

• Thin-layer viscous in 3 directions - two wall corner model 

- CFL3DE 

• Streamwise-relaxation crossflow-AF, space-marching Euler 
or PNS, first or second order 

® Perfect gas or equilibrium air 


INDUSTRY USE STATUS OF CODES 


NASP prime contractors - all instructed in code use at workshop 
August 2-3, 1988. 

- Production code for NASP for McDonnell Aircraft Company 

- Some use at General Dynamics 
Other industry use 

- Boeing 

- Northrup 

- United Technologies Research Center 
University use 

- Iowa State University 
Other government use 


Naval Surface Warfare Center 


SR71 - GEOMETRY DEFINITION 
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SR71 - EQUILIBRIUM AIR EULER SOLUTION 


CROSS-SECTION PRESSURE DISTRIBUTIONS AND GRID DEFINITIONS 
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GRID DENSITY COMPARISON 


CROSS-SECTION 45 (COWL PLANE) 



65 X 65 X 46 GRID 
Y-PLUS ORDER 1. 


SOX 130X56 GRID 
Y-PLUS ORDER .001 



PRESSURE COMPARISONS 


M = 12.55 

Re = 2.7 million/ft. 

Zero degrees angle of attack 



Fuselage Station (inches) Fuselage Station (inches) 

Open symbols - upper surface 
Closed symbols - lower surface 


HEAT TRANSFER COMPARISONS 


M = 12.55 

Re = 2.7 million/ft. 

Zero degrees angle of attack 



Fuselage Station (inches) 



Open symbols - upper surface 
Closed symbols - lower surface 


AERODYNAMIC COEFFICIENTS COMPARISON MCDONNELL 
BLENDED WING BODY CFD AND EXPERIMENT 

M = 12.4, a = 6°, Re/L = 930,000/ft 
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DENSITY 

GENERIC NOZZLE AFTERBODY 
^rCONPUTATIONAL METHODS BRANCH NASA LANGLY 
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FUTURE PLANS RELATED TO THE NATIONAL AERO-SPACE 

PLANE PROGRAM 


Continue to expand the envelope of capabilities for the code to include 
calculations of an entire NASP-like configuration 


- Improved zonal capabilities for inlets with sweep and combustors 


- Addition of non-equilibrium chemistry for combustor and 
nozzle/afterbody calculations 


- All capabilities scheduled for production code by 1/90 (NASP 
Technology Maturation Program) 


SESSION VIII 


SPACE SHUTTLE 


Chairman: 

Luke A. Schutzenhofer 

Chief, Computational Fluid Dynamics Branch 
NASA Marshall Space Flight Center 
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NUMERICAL AERODYNAMIC SIMULATION OF THE 
SPACE SHUTTLE ASCENT ENVIRONMENT 


J.P. Slotnick 

Lockheed Engineering and Sciences Company, Houston, TX 


F.W. Martin. Jr. 

NASA Lyndon B. Johnson Space Center, Houston, TX 


in collaboration with the 


NASA Ames Space Shuttle Flow Simulation Group 
P.G. Buning, I.T. Chiu. R.L. Meakin, S. Obayashi. Y.M. Rizk 
S. Ben-Shmuel. J.L. Steger. and M. Yarrow 
NASA Ames Research Center, Moffett Field, CA 


ABSTRACT 


After the STS 51-L accident, an extensive review of the Space Shuttle Orbiter’s 
ascent aerodynamic loads uncovered several questionable areas that required fur- 
ther analysis. The insight gained by comparing the Shuttle ascent CFD numerical 
simulations, obtained by the NASA Ames Space Shuttle Flow Simulation Group 
led by Dr. J.L. Steger, to the current IVBC-3 aerodynamic loads database was 
instrumental in resolving uncertainties on the Orbiter payload bay doors and fuse- 
lage. Initial confidence in the numerical simulations was gained by comparing them 
with the limited flight data that had been obtained during the Orbiter Flight Test 
(OFT) program. Current CFD results exist for Mach numbers 0.6, 0.9, 1.05, 1.55, 
2.0, and 2.5. Since the pre STS-1 wind tunnel test program (IA-105) often yields 
considerable differences when compared to STS-5 flight data, the Moo=1.05 tran- 
sonic case is the most investigated. The IA308 mated-vehicle hot gas plume wind 
tunnel test, recently completed at AEDC 16T (transonic) and Lewis (hypersonic), 
is also used to compare with the computation where applicable. 
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Objective 


o Joint effort between the ARC team led by J.L. Steger and P.G. Bun- 
ing and a JSC team led by F.W. Martin, Jr. 

o ARC Team : Develop the technology to numerically simulate 

the complex launch vehicle (SSLV) geometry in 
a time accurate manner to 

• compute SSLV flowfield data 

• support Fast Separation Abort study 

• assist JSC with Orbiter problematic issues: 

• Crew Escape flowfield 

• STS-27 Debris Study flowfields 

o JSC Team : Apply the CFD technology to gain insight into 

the ascent aerodynamic loads environment 

• SSLV payload bay pressure distribution 

• Orbiter Wing Loads 


Jeff P. Slot nick / March 8, 1989 


Scope 


o ARC Team 

• Model the Orbiter ascent configuration with the Chimera 
composite grid discretization approach. 

• Overset body-conforming grids for each major component 
(CALS PAN Pegasus code). 

• Apply F3D , an implicit approximately factored finite- 
difference procedure, to solve the three-dimensional thin- 
layer Navier-Stokes equations. 

• Employ the NAS Cray 2 and Cray Y-MP. 


o JSC Team 

• Provide technical guidance and evaluate ARC results with 
respect to existing databases: 

t> I A- 105 WTT, high resolution model. 

> IA-308 WTT, hot gas plume simulation. 

> Limited flight data from STS-1 through STS-5. 

> IVBC-3 operational database. 


Jeff P. Slotnick / March 8, 1989 


Justification 


The Orbiter’s modest structural capability, coupled with a complex launch 
vehicle configuration and severe ascent environment, places unprecendented 
demands on the aerodynamacist. 

o Orbiter structure was designed to civil transport specifications. 

• Consistent with entry flight envelope. 

• q « 350 psf, g 2 =2.5 (g z =2.25 as weight increases). 

o Ascent environment is considerably more severe, q=750-819 psf, 
with the simultaneous occurance of: 

• maximum q 

• jet stream penetration (perturbates a,/3 and q) 

• transonic flow field 

• maximum differential pressure, inside-outside, as the vehi- 
cle vents 

• plume/flow field interactions 

• High R e 100 million) at maximum load condition 


Jeff P. Slotnick / March 8. 1989 
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Historical Perspective 


o Discrepancies exist between aerodynamic predictions and flight ex- 
perience. 



4 • 

MACH NUMBER 



PITCHING MOMENT CAUSES THE SSV TO 
FLY AT A STEEPER ATTITUDE RESULTING 
IN AN INCREASE IN THE VERTICAL 
COMPONENT OF THRUST 


NORMAL FORCE CAUSES AN INCREASE IN 
NORMAL ACCELERATION A THE LOAO 
RELIEF SYSTEM COMMANOS A NOSE 
DOWN PITCH 


o Force and moment data was easily corrected with flight derived 
aerodynamic increments. 

o Aerodynamic loads (pressure distribution) cannot be readily cor- 
rected because of limi ted flight pressure measurements. 
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Current Database 


o The /ntegrated Vehicle Baseline Configuration (IVBC) - 3 aero loads 
database is a composite of 

• IA-105 WTT data 

• Flight pressure measurements 

• Flight strain measurements 

• Engineering judgement 


o This resulting database contains relatively large uncertainties which 
have led to structural analysis predictions that were inconsistent 
with flight experience, especially in the payload bay door area. 


Jeff P. Slotnick / March 8, 1989 


Selected STS-5 Tap Points 



Jeff P. Slot nick / March a, 1989 


124 


AERO DATABASE COMPARISONS 

TAP V07P9108A CP HISTORY 



AERO DATABASE COMPARISONS 

TAP V07P9495A CP HISTORY 
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AERO DATABASE COMPARISONS 

ORBITER Upper Pitch Plane CP Distributions 



AERO DATABASE COMPARISONS 

ORBITER Mid-Fuselage Region CP Distributions 
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PAYLOAD BAY DOOR PRESSURE COMPARISONS 
ORBITER X-STATION (XO) = 590 inches 
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Status 


o Preliminary CFD results have been instrumental in gaining insight 
into the present design aero database. Even with the large ad- 
vances in SSLV geometry modelling, all three aero data sources 
(flight, wind tunnel, and CFD) are required to completely under- 
stand Shuttle ascent aerodynamics and structural loads. 

o On-going CFD Activities: 

• Add the SRB attach ring and IEA box protuberances into 
the SSLV grid system. 

• Add more realistic aft attach hardware to the geometry 
model to better simulate blockage effects. 

• Add 7=7(T) to the F3D flow solver. 

• Grid the vertical tail and include in the geometry model. 

• Grid the SSMEs and incorporate the SSME plumes into 
the solution. 
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COMPUTATIONAL FLUID DYNAMICS ANALYSIS OF SPACE SHUTTLE MAIN 
PROPULSION FEED LINE 17 -INCH DISCONNECT VALVES 

Max Kandula* and Daniel Pearce** 

Lockheed Engineering and Sciences Company 
Houston , Texas 


ABSTRACT 


A steady incompressible three-dimensional (3-D) viscous flow 
analysis has been conducted for the Space Shuttle Main 
Propulsion External Tank (ET)/Orbiter (ORB) propellant feed line 
quick separable 17-inch disconnect flapper valves for liquid 
oxygen (L0 2 ) and liquid hydrogen (LH 2 ) . The main objectives of 
the analyis were to predict and correlate the hydrodynamic 
stability of the flappers and pressure drop with available water 
test data. 

Computational Fluid Dynamics (CFD) computer codes were procured 
at no cost from the public domain, and were modified and extended 
to carry out the disconnect flow analysis. The grid generator 
codes SVTGD3D and INGRID, developed by Sverdrup Technology Inc. , 
were obtained from Arnold Air Force Station, Tennessee. NASA 
Ames Research Center supplied the flow solution code INS3D, and 
the color graphics code PL0T3D. A driver routine was developed 
to automate the grid generation process. Components such as 
pipes, elbows and flappers can be generated with simple commands, 
and flapper angles can be varied easily. The flow solver INS3D 
code was modified to treat interior flappers, and other 
interfacing routines were developed, which include a turbulence 
model, a force/moment routine, a time-step routine, and initial 
and boundary conditions. In particular, an under-relaxation 
scheme was implemented to enhance the solution stability. 

Major physical assumptions and simplifications made in the 
analysis include the neglect of linkages (drive/follower arms) , 
slightly reduced flapper diameter, and smooth solid surfaces. A 
grid size of 54x21x25 was employed for both the L0 2 and LH 2 
units. Mixing length theory applied to turbulent shear flow in 
pipes formed the basis for the simple turbulence model. 

Results of the analysis are presented for L0 2 and LH 2 
disconnects. The predicted stop loads, hydrodynamic stability 
boundaries of the ET and orbiter flappers, and pressure drop 
across the valve compare well with the water test data, covering 
a tube Reynolds number of 3.5E06 for L0 2 unit and 2.4E06 for LH 2 
unit. The ability to predict the valve performance and flapper 
stability over a wide range of flow rates and flapper angle 
combinations demonstrates the validity of the CFD model. 


* Advanced Systems Engineering Specialist 

** Engineer 
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lh 2 disconnect valve 



Re = 2.7E06 50E06 


OBJECTIVES 


PREDICT 

• HYDRODYNAMIC STABILITY 
. PRESSURE DROP 




CAVITATION 
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COMPUTATIONAL GRID 


CODE : SVTGD3D ( SONI ET AL. 1986 ) 
ORG : SVERDRUP / AIR FORCE 



L02 DISCONNECT 


COMPUTATIONAL GRID 

L02 DISCONNECT 


FLOW SOLUTION DEVELOPMENT 


CODE : INS3D ( KWAK ET AL. , 1985 ) 
ORG : NASA AMES 


• INTERNAL OBSTACLES 

• TURBULENCE MODEL 
. UNDER-RELAXATION 

ANGLE CORRECTIONS 



f ET : FA = (PA+1) deg 
) . J ORB : FA = (PA-5) deg 

2 ] PA = PRELOAD ANGLE 


FA = FLIGHT ANGLE 


LH •• NOME 
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UNDER-RELAXATION 


• ENHANCED STABILITY 

n+1 n n 

Qj = Qj + ccr AQj 

Qi = [p,u,v,w] T 

ccr = RELAXATION FACTOR 
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EFFECT OF FLOW RATE 
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STABILITY BOUNDARIES 
L0 2 DISCONNECT 
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STABILITY BOUNDARIES 
LH 2 DISCONNECT 
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SUMMARY 


. ENHACED THE STABILITY OF INS3D 

. VALIDATED THE CODE WITH DATA 

. INTERNAL FLOW 

• INTERNAL OBSTACLES 

• HIGH REYNOLDS NUMBER (3.5E06) 

. PREDICTED FLAPPER STABILITY 

. ASSISTED IN THE VALVE REDESIGN 

• LATCH BETWEEN FLAPPERS 
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Analysis of the SSME HPOTP Bearing Inlet Cavity 


P. McConnaughey , NASA/MSFC 


Analysis of the flow in the Space Shuttle Main Engine (SSME) 
high pressure oxygen turbopump (HPOTP) bearing # 1 inlet cavity 
has been completed in support of return-to-f light . With the 
incorporation of several design changes in the Phase II 
turbopump, rotordynamic stability of the pumps has been 
enhanced, but the durability and life of the LOX-cooled 
bearings has decreased. During the post-Challenger SSME 
recertification, the MSFC bearing team investigated the causes 
of limited bearing durability. One topic addressed by this 
team was the flow environment upstream of the pump-end bearing 
and the effect of seal exit swirl (Phase I labyrinth seal vs. 
Phase II damping seal) and a cavity anti-vortex rib on the 
bearing environment and life. The objective of the present 
work was to define the hydrodynamic environment upstream of the 
pump-end bearing and determine the effect of seal exit swirl 
and the anti-vortex rib on bearing inlet swirl. 

The problem was posed as an axisymmetric cavity flow with the 
computational domain extending from the seal exit to the 
bearing inlet. This domain was discretized with 22800 grid 
points. Boundary conditions were obtained from a 1-D model of 
the SSME coolant path. These resulted in an axial Reynolds 
number of 297000 and with a seal tip speed of 29,200 rpm. The 
inlet Mach number was 0.19 and the problem was solved with the 
CMINT code utilizing the Briley-McDonald/Beam-Warming algorithm 
with preconditioning to speed convergence at low Mach numbers. 
Three parametric cases with inlet swirl of 50% shaft speed 
(labyrinth seal), 20% shaft speed (damping seal) and no inlet 
swirl were considered. Total CPU time for all analyses was 9 
hours on a Cray X-MP with memory requirements of 1.7 million 
words . 

Computational results indicate large vortical flow structures 
in the cavity, with the labyrinth, damping, and no-swirl cases 
yielding bearing inlet swirl rates of 14, 10, and 9 percent of 
shaft speed, respectively. These small differences are due to 
fluid spin-up on the shaft and inner race and indicate that 
upstream influences (either inlet swirl or anti-vortex ribs) 
have little effect on bearing inlet conditions. When these 
results were used as input to the SHABRETH bearing model, 
limited durability could not be explained by these small 
differences in swirl. Also, based on these results, a proposed 
design change for the cavity anti-vortex rib was not 
implemented by the SSME chief engineer. 
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BEARING INLET CAVITY ANALYSIS 


OVERVIEW 

• OBJECTIVE 

• JUSTIFICATION 

• APPROACH 

• RESULTS 

• SUMMARY AND CONCLUSIONS 

• PROGRAM IMPACT 

OBJECTIVES (IN SUPPORT OF RETURN-TO-FLIGHT) 

• DEFINE THE HYDRODYNAMIC ENVIRONMENT FOR THE LOX 
BEARING #1 INLET CAVITY 

- UNDERSTAND THE SEAL EXIT/CAVITY/BEARING INLET FLOW PATH 

- QUANTIFY THE PRESSURE DROP IN THE CAVITY 

• EVALUATE THE EFFECT OF SEAL EXIT SWIRL ON BEARING INLET SWIRL 

- INLET SWIRL AFFECTS VISCOUS HEAT GENERATION IN BEARING, CAGE 
DRAG AND PRESSURE DROP ACROSS BEARINGS 

• EVALUATE EFFICACY OF THE BEARING CAVITY ANTI-VORTEX RIB 
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BEARING INLET CAVITY ANALYSIS 


JUSTIFICATION 

• DURABILITY OF HPOTP PUMP-END BEARINGS (2000 SECOND LIMIT) 

- LOWER LIFE IN PHASE II DESIGN RELATIVE TO PHASE I DESIGN 

- CONFIGURATION CHANGES 

• DAMPING SEALS (PHASE II) VS. LABYRINTH SEALS (PHASE I) 

• INCREASED DEADBAND, INTERNAL CLEARENCES 

• ELONGATED CAGE POCKETS 

•PROPOSED DESIGN CHANGE TO REMOVE ANTI-VORTEX RIB IN CAVITY 
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BEARING INLET CAVITY ANALYSIS 


PHASE I 


PHASE II 


ANTI-VORTEX 

RIB 
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APPROACH 


BEARING INLET CAVITY ANALYSIS 


PROBLEM DEFINITION 

• DOMAIN FROM SEAL EXIT TO BEARING INLET 

- AXISYMMETRIC ANALYSIS (NO ANTI-VORTEX RIB) 

- LEAKAGE PATH IGNORED (1% OF FLOW) 

- ALL FLOW TO BEARING GOES BETWEEN INNER RACE AND CAGE 

• Re ^ = 297000, Q = 29200 rpm, M j n = 0.1 9 (1 09% POWER LEVEL) 

METHOD OF ANALYSIS 

• COMPRESSIBLE NAVIER-STOKES EQUATIONS (2-D AXISYMMETRIC WITH SWIRL) 

• BRILEY-McDONALD/BEAM-WARMING WITH PRECONDITIONING 

• CM I NT CODE 

• LOW, Re K- e TURBULENCE MODEL 

• 22800 GRID POINTS 

• SOFT-START ON DOWNSTREAM STATIC PRESSURE 

• 9 HOURS, 1 .7 X 1 0 6 WORDS ON MSFC CRAY X-MP 

VALIDATION OF DISK CAVITY ANALYSES 

• OWEN AND PINCOMBE (JFM, 1980) 


• DAILY AND NECE (1960) 


SSME HPQTP BEARING #1 INLET CAUITY ANA1 YSIS 
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LOX SWIRL PROFILES AT BEARING CAGE INLET 
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BEARING INLET CAVITY ANALYSIS 
SUMMARY AND CONCLUSIONS 


• HYDRODYNAMIC ENVIRONMENT FOR THE SSME HPOTP BEARING #1 
INLET CAVITY HAS BEEN DEFINED FOR PHASE I AND PHASE II DESIGNS 

• SMALL DIFFERENCE (4%) IN BEARING INLET SWIRL BETWEEN PHASE I 
AND PHASE II PUMPS 

Ln 

oo 

• CFD RESULTS WERE USED AS INPUT INTO SHABRETH BEARING MODEL 
IN SUPPORT OF RETURN-TO-FLIGHT 
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BEARING INLET CAVITY ANALYSIS 


PROGRAM IMPACT 


• DETERMINED THAT THE LIMITED DURABILITY OF THE HPOTP PUMP-END 
BEARINGS IS NOT CAUSED BY CHANGES IN SWIRL ASSOCIATED WITH THE 
PHASE II DESIGN DAMPING SEALS 

• DETERMINED THAT THE PROPOSED MODIFICATION OF THE CAVITY ANTI- 
VORTEX RIB WOULD HAVE LITTLE EFFECT ON SWIRL ENTERING THE BEARING. 
PROPOSED DESIGN CHANGE WAS NOT IMPLEMENTED BY CHIEF ENGINEERS 
OFFICE 
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o DISK CAVITY ANALYSIS AT NASA-MSFC 
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A COMBINED EULERIAN-LAGRANGIAN TWO-PHASE 
ANALYSIS OF THE SSME HPOTP NOZZLE 
PLUG TRAJECTORIES 


R. Garcia, P. K. McConnaughey , NASA/MSFC 
F. de Jong, J. Sabnis, Scientific Research Associates 
D. Pribik, Rocketdyne Division, Rockwell International 


As a result of high cycle fatigue, hydrogen embrittlement, 
and extended engine use, it was observed in testing that the 
trailing edge on the first stage nozzle plug in the High 
Pressure Oxygen Turbopump (HPOTP) could detach. The 
objective of this study was to predict the trajectories 
followed by particles exiting the turbine. Experiments had 
shown that the heat exchanger coils, which lie downstream of 
the turbine, would be ruptured by particles traveling in the 
order of 360 ft/sec. An axisymmetric solution of the flow 
was obtained from the work of Lin et al [1] who used INS3D to 
obtain his solution. The particle trajectories were obtained 
using the method of de Jong et al [2] which employs 
Lagrangian tracking of the particle through the Eulerian flow 
field. The collision parameters were obtained from 
experiments conducted by Rocketdyne using problem specific 
alloys, speeds, and projectile geometries. A complete 3-D 
analysis using the most likely collision parameters shows 
maximum particle velocities of 200 ft/sec. in the heat 
exchanger region. Subsequent to this analysis, an engine 
level test was conducted in which seven particles passed 
through the turbine but no damage was observed on the heat 
exchanger coils. 


1) Lin, S.-J. and Chang, J.L.C.: "Numerical Study of Laminar 
and Turbulent Flow Inside a Turnaround Duct With and Without 
Guide Vanes," AIAA Paper 87-0365, 1987. 

2) de Jong, F.J., Sabnis, J.S., and McConnaughey, P.K.: M A 
Combined Eulerian-Lagrangian Two-Phase Flow Analysis of SSME 
HPOTP Nozzle Plug Trajectories; Part I- Methodology," To be 
presented at the AIAA Joint Propulsion Conference, 1989. 
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HIERARCHY OF ANALYSIS 


® 2-D: PARTICLE EXITS TURBINE WITH NO SWIRL 

• 3-D: PARTICLE WITH SWIRL, NO RADIAL STRUTS 

• 3-D: PARTICLE WITH SWIRL, RADIAL STRUTS INCLUDED 


PARAMETERS OF ANALYSIS 

• 15 DIFFERENT PARTICLE RADIAL POSITIONS AT THE INLET TO THE DOMAIN FOR EACH 
RUN. (INNER TO OUTER WALL) 

• DRAG COEFFICIENT OF PARTICLES - 0.4, 0.6, 0.8 

• TURBINE EXIT AXIAL VELOCITY OF PARTICLE FROM 180 TO 255 FT/SEC. 

• SWIRL VELOCITY OF PARTICLE VARIED FROM ZERO TO BLADE SPEED 

• LOSS OF NORMAL MOMENTUM AT THE WALLS: 

• COEFFICIENT OF RESTITUTION - 0.2, 0.4, 0.6, 0.8 

• LOSS OF TANGENTIAL MOMENTUM AT THE WALLS: 

• COEFFICIENT OF SLIDING FRICTION - 0.0, 0.1 , 0.2, 0.3 
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SSME HPOTP NOZZLE PLUG TRAJECTORIES 
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OBJECTIVE 


• TO DETERMINE THE LOCATIONS AND VELOCITIES OF NOZZLE PLUG IMPACTS ON 
THE HPOTP LOX HEAT EXCHANGER. TEST DATA INDICATED IMPACT VELOCITIES OF 
360 FT/SEC WOULD CAUSE HEAT EXCHANGER COIL RUPTURE. 


JUSTIFICATION 

• RESULTS REQUESTED BY SSME CHIEF ENGINEER TO ASSESS: 

• RISK OF FULL SCALE ENGINE TEST. 

• RISK OF RETURN TO FLIGHT OF THE SPACE SHUTTLE WITH CURRENT DESIGN. 
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TECHNICAL APPROACH 


• USED UN’s AND CHANG’S (ROCKETDYNE) AXISYMMETRIC SOLUTION OF THE FLOW 
IN THE HPOTP HOT GAS MANIFOLD 

• SOLUTION OBTAINED USING INS3D: 

• INCOMPRESSIBLE M = 0.1 

• Re =3 x 1 0 6 , 588 FT/SEC AXIAL VELOCITY 

• TRAJECTORIES CALCULATED USING LAGRANGIAN TRACKING OF THE PARTICLES 
THROUGH THE EULERIAN FLOW FIELD 

• METHOD OF F. J. de JONG AND J. S. SABNIS OF SRA 

• F = ma ON THE PARTICLE THROUGH THE FLOW FIELD 

• COLLISION MODEL AT THE WALLS. 

• VALIDATION 

• INS3D: COMPARED TO EXPERIMENTAL MEASUREMENTS FOR 
VARIOUS INCOMPRESSIBLE INTERNAL FLOWS. 

• LAGRANGIAN TRACKER: COMPARED AGAINST EXPERIMENTAL 
MEASUREMENTS IN A TWO-PHASE MIXING LAYER FLOW. 
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SSME HPOTP NOZZLE PLUG TRAJECTORIES 


OVERVIEW 

• INTRODUCTION TO THE PROBLEM 

• OBJECTIVES 

• TECHNICAL APPROACH 

• RESULTS 

• CONCLUSION/IMPACT 
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3-D ANALYSIS WITH RADIAL STRUTS (COMPLETE GEOMETRY) 
PARTICLE VELOCITIES* AT HEAT EXCHANGER (FT/SEC) 
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IMPACT 


• DECISION MADE TO CONDUCT FULL SCALE TEST WITH PURPOSELY 
CUT NOZZLE PLUGS: 

NO DAMAGE TO THE HEAT EXCHANGER 

• DECISION MADE TO USE CURRENT DESIGN ON RETURN TO FLIGHT 
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CONCLUSIONS 


• PARTICLE MOTION DOMINATED BY COLLISION PARAMETERS 

• MOST PARTICLES WILL HIT THE HEAT EXCHANGER DUE TO “FOCUSING” OF THE TURNING 
VANES. 

• MEAN COLLISION PARAMETERS (C RES = 0.6, n = 0.1) YIELD MAXIMUM PARTICLE IMPACT 
VELOCITIES OF 180 FT/S 

• WORST CASE STACK OF COLLISION PARAMETERS YIELDS MAXIMUM IMPACT VELOCITIES 
OF 300 FT/S (RUPTURE VELOCITY IS 360 FT/S). 
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Conjugate ( Solid/Pluid ) Computational Fluid Dynamics 
Analysis of the Space Shuttle Solid Rocket Motor 
Nozzle/Case and Case Field Joints 


D. Doran, NASA/MSFC, L. W. Keeton, P. J. Dionne, 
and A. K. Singhal, CFD Research Corporation 


This work describes three-dimensional, conjugate 
(solid/fluid) heat transfer analyses of new designs of the 
Solid Rocket Motor (SRM) nozzle/case and case field joints. 

The main focus of the study has been to predict the 
consequences of multiple "rips" (or debonds) in the ambient 
cure adhesive packed between the nozzle/case joint surfaces 
and the bond line between the mating field joint surfaces.' 
The models calculate the transient temperature responses of 
the various materials neighboring postulated f low/leakpaths 
into, past and out from the nozzle/case primary 0-ring cavity 
and case field capture 0-ring cavity. These results were 
used to assess if the design was failsafe (i.e. no potential 
0-ring erosion) and reusable (i.e. no excessive steel 
temperatures ) . 

The models are adaptions and extensions of the general 
purpose PHOENICS fluid dynamics code. 

A non-orthogonal coordinate system was employed and 11,592 
control cells for the nozzle/case and 20,088 for the case 
field joints are used with non-uniform distribution. 

Physical properties of both fluid and solids are temperature 
dependent . 

A number of parametric studies were run for both joints with 
results showing temperature limits for reuse for the steel 
case on the nozzle joint being exceeded while the steel case 
temperatures for the field joint were not. 0-ring 
temperatures for the nozzle joint predicted erosion while for 
the field joint they did not. 
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SRB NOZZLE/CASE AND CASE FIELD JOINT 
FLOW AND THERMAL ANALYSIS 


OVERVIEW 


| 

• PROBLEM DEFINITION 

• OBJECTIVES 

• TECHNICAL APPROACH 

• RESULTS 

• PROGRAM IMPACT AND CONCLUSIONS 
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SRB NOZZLE/CASE AND CASE FIELD JOINT 
FLOW AND THERMAL ANALYSIS 


OBJECTIVE: 


JUSTIFICATION: 


DEFINITION OF THE FLOW AND THERMAL ENVIRONMENTS IN THE 
SRB NOZZLE/CASE AND CASE FIELD JOINTS ASSUMING SEVERAL 
POSTULATED DOUBLE LEAKPATHS THROUGH THE NOZZLE/CASE 
AMBIENT CURE ADHESIVE AND THE CASE FIELD BOND LINE 
BETWEEN THE MATING FIELD JOINT SURFACES 


TO DEFINE CONSERVATIVE, BUT CREDIBLE DESIGN AND FAIL-SAFE 
CONDITIONS AND CRITERIA FOR NOZZLE/CASE AND CASE FIELD 
JOINT THERMAL ANALYSIS 


DESIGN - TEMPERATURES WHICH GUARANTEE REUSABILITY 

OF METAL PARTS AND NO EROSION OF O-RING SEALS 


FAIL-SAFE - TEMPERATURES WHICH GUARANTEE A STRUCTURAL 
SAFETY FACTOR OF 1 .0 ON METAL PARTS AND NO 
EROSION OF ONE SEAL 
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SRB NOZZLE/CASE AND CASE FIELD JOINT 
FLOW AND THERMAL ANALYSIS 


APPROACH: APPLICATION OF 3-D NAVIER-STOKES CODE (PHOENICS) TO FLOW 

IN THE DEBOND GAP AND O-RING REGIONS COUPLED WITH 
HEAT CONDUCTION OF THE FLUID AND NEIGHBORING “SOLIDS” 
(NOZZLE/CASE 1 1 592 CELLS; CASE FIELD 20088 CELLS) 


PHOENICS 

00 

• (PARABOLIC, HYPERBOLIC, OR ELLIPTIC NUMERICAL 
INTEGRATION CODE SERIES) 

MODEL USES: 

• BODY-FITTED COORDINATES 

• FINITE-VOLUME FORMULATION 

• 3-D TRANSIENT 

• K-EPSILON TURBULENCE MODEL 

• TEMPERATURE DEPENDENT SOLID AND FLUID 
PROPERTIES 

• WALL FUNCTIONS (LOG LAW OF THE WALL) 

• CHILTON-COLBURN FORM OF REYNOLDS 
ANALOGY FOR HEAT TRANSFER AT WALL 
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SUMMARY OF MAXIMUM TEMPERATURES FOR 
NOZZLE/CASE JOINT MODEL 


Run 

No. 

Insul. 

Gap 

(in.) 

Steel 

Gap 

(in.) 

Debond 
Width (5) 
(in.) 

Diff. P 
(psi) 

Diff. P 
Angle 
(deg) 

T max 
Gas 
(deg F) 

T max 
Steel 
(deg F) 

T max 
O-ring 
(deg F) 

Notes 

1 

.024 

.004 

F.O. 

1.08 

0-180 

142 

180 

105 


2 

.024 

.010 

F.O. 

1.08 

0-180 

295 

230 

170 


3 

.024 

.004 

3.5 

1.08 

0-180 

1120 

1092 

769 


4 

.024 

.010 

3.5 

1.08 

0-180 

2071 

1432 

1423 


5 

.024 

.015 

F.O. 

envelope 

0-180 

258 

132 

129 


6 

.100 

.050 

F.O. 

envelope 

0-180 

1125 

659 

607 


7 

.024 

.008 

3.5 

envelope 

0-180 

1628 

1117 

885 


8 

.024 

.008 

3.5 

QM-4 

0-180 

1103 

742 

582 


9 

.024 

.008 

3.5 

QM-4 

0-180 

1153 

769 

601 

1 

10 

.024 

.008 

3.5 

QM-4 

0-120 

1256 

852 

667 

1 

11 

.024 

.008 

3.5 

QM-4 

0-120 

1643 

1380 

909 

1,2 

12 

.050 

.010 

0.4 

QM-4 

0-120 

2898 

2320 

1510 

1,2 

13 

.050 

.010 

0.4 

QM-4 

0-120 




1,2,3 

14 

.024 

.008 

3.5 

QM-4 

0-120 

1502 

1304 


1,2,3 

15 

.046 

.008 

0.75 

QM-4 

0-120 

2276 

2352 

1856 

1,2,3 

16 

.090 

.008 

0.75 

QM-4 

0-120 

2834 

2846 

2341 

1,2,3 

17 

.046 

.008 

0.75 

QM-4 

0-120 

2278 

2359 

1859 

1,2, 3, 4 

18 

.046 

.008 

0.75 

QM-4 

0-120 

2278 

2351 

1858 

1,2, 3, 6 

19 

.046 

.008 

0.75 

QM-4 

0-120 

2211 

2252 

1847 

1,2, 3, 7 

20 

.046 

.008 

0.75 

QM-4 

0-120 

2083 

2353 

1691 

1,2, 3, 8 

21 

.046/. 020 

.008 

0.75 

QM-4 

0-120 

2089 

2178 

1588 

1,2,3, 9 

22 

.046/.020 

.008 

0.75 

QM-4 

0-120 

2059 

2187 

1568 

1,2,3,9,10 

23 

.046/.020 

.008 

0.75 

QM-4 

0-120 

2053 

2272 

1543 

1,2, 3, 9, 11 

24 

.046/.020 

.008 

0.75 

QM-4 

0-120 

1866 

2255 

1509 

1,2,3,9,11,12,13 

25 

.046/. 020 

.008 

0.75 

DM-8 

0-120 

1610 

1818 

1274 

2,3,9, 11 


Notes: 1 . Improved modeling of QM-4 duty cycle 

2. Model revised to open NBR/glass phenolic area 360 deg. 
and improve average material properties 

3. Material properties computed as function of temperature 

4. Same as Run 15 but with number of time steps doubled 

5. One-half of values listed input to model due to symmetry 

6. 2 X nominal number of grid nodes in Z direction 


7. Increased number of grid nodes in X direction (<2 X nominal) 

8. 2 X nominal number of grid nodes in X direction 

9. Model regrided to provide two insulation gap width capability 

10. Y grid sensitivity. Number of insulation cells increased 

11. K-Epsilon turbulence model 

12. Y grid sensitivity. Number of fluid cells increased by factor of 2 

13. Calculation of O-ring wall temperature corrected 
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NOZZLE/CASE JOINT PHOENICS RUN 
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SUMMARY OF MAXIMUM TEMPERATURES 
FOR FIELD JOINT MODEL 


RUN 

SPECIFICATION 

CALCULATED RESULTS 

NO. 

G1 

(IN) 

CIRCUM. 

DEBOND 

WIDTH 

(IN) 

CIRCUM. 

DEBOND 

ANGLE 

(DEG) 

Ap 

PSI 

TMAX 
STEEL 
(DEG C) 

TMAX 
GAS 
(DEG C) 

TMAX 
O-RING 
(DEG C) 

1 

0.002 

2.8 

120 

0.5 

16.1 

13.9 

13.6 

2 

0.002 

2.8 

15 

0.5 

16.2 

14.1 

13.7 

3 

0.016 

0.25 

120 

0.5 

16.5 

23.3 

14.2 

4 

0.016 

0.25 

15 

0.5 

161.7 

746.6 

607.4 

5 

0.016 

0.25 

15 

0.1 

21.0 

65.4 

42.3 


NOTE: FOR ALL CASES G2 = 0.0264", AND G3 - 0.0241" 
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SRB NOZZLE/CASE AND CASE FIELD JOINT 
FLOW AND THERMAL ANALYSIS 


PROGRAM IMPACT: 


CONCLUSIONS: 


• STEEL CASE TEMPERATURE LIMIT FOR REUSE (1000°F) 

• EXCEEDED FOR NOZZLE/CASE JOINT (IF 2 LEAKPATHS OCCUR WHICH IS 
VERY UNLIKELY) 

• WITHIN LIMIT FOR CASE FIELD JOINT 

• ABLATION TEMPERATURE LIMIT FOR O-RINGS (700°F) 

• EXCEEDED FOR NOZZLE/CASE JOINT (IF 2 LEAKPATHS OCCUR WHICH 
IS AGAIN UNLIKELY) 

• WITHIN LIMIT FOR CASE FIELD JOINT (EXCEPT ONE CASE WHERE 
FLOW PATH IS ONLY OVER 15°) 


• ANALYSIS CONSIDERED CONSERVATIVE, DUE TO WORST CASE 
SCENARIOS BEING ANALYZED 

• ANALYSIS PERFORMED TO VERIFY OTHER ANALYSIS BEING DONE IN 
PARALLEL, AT MORTON THIOKOL TO SUPPORT RETURN-TO-FLIGHT 
STATUS 

• PROBABILITY OF TWO LEAKPATHS OCCURING IS ONE IN ONE-THOUSAND 
SO DESIGN MEETS FAIL-SAFE AND REUSEABILITY CONDITIONS. 
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SIMULATION OF TURBOMACHINERY FLOWS 
by 

John J. Adamczyk 
NASA Lewis Research Center 

ABSTRACT 

Significant advancements have been made in the last five years in the abil- 
ity to model turbomachinery flows of engineering interest. This advance- 
ment can be directly attributed to the second generation of supercomputers 
like the Cray XMP and Cray II and advanced instrumentation techniques. 

Early on, the National Aeronautics and Space Administration Lewis Research 
Center recognized the potential gains in turbomachinery performance and 
life that could be achieved by taking advantage of this technology and 
instituted a comprehensive research program in turbomachinery flow model- 
ing. This activity combined the areas of fluid flow analysis, computa- 
tional fluid dynamics, and experimental fluid mechanics. As a result of 
this activity, Lewis has become an internationally recognized leader in 
turbomachinery flow modeling. Many of the research activities conducted 
under this progam have been utilized by industry. The presentation will 
give an overview of this program and provide sample illustration of simu- 
lation performed to date. 
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PREDICTION OF TURBINE ROTOR-STATOR INTERACTION 
USING NAVIER-STOKES METHODS 

N. K. Madavan and M. M. Rai 
NASA Ames Research Center 
Moffett Field, CA 94035 

and 

S. Gavali 

Amdahl Corporation 
Sunnyvale, CA 94088 


Abstract 


Flows in turbomachinery are generally complex and do not easily lend themselves to 
numerical computation. The flows are three-dimensional and inherently unsteady. The 
unsteadiness arises from the interaction of the downstream airfoils with the wakes and 
passage vortices generated upstream, from the motion of the rotors relative to the stators 
and from vortex shedding at blunt airfoil trailing edges. Complicated blade geometries and 
flow phenomena such as separation and periodic transition from laminar to turbulent flow 
add to the numerical complexity. Nevertheless, the accurate numerical analysis of such 
flows is a problem of considerable interest and practical importance to the turbomachinery 
community. 

Much of the early work in turbomachinery flow prediction focussed on airfoil cascades. 
While such analyses of flows in isolated airfoil rows have helped improve our understanding 
of the flow phenomena and have gained widespread acceptance in the industrial community 
as a design tool, they do not yield any information regarding the unsteady effects arising 
out of rotor-stator aerodynamic interaction. These interaction effects become increasingly 
important as the distance between successive stator and rotor rows is decreased. Thus, 
the need exists for analytical tools that treat the rotor and stator airfoils as a system and 
provide information regarding the magnitude and the impact of the unsteady effects. 

The focus of this presentation is a three-dimensional, time- accurate, thin-layer Navier- 
Stokes code that has been recently developed to study rotor-stator interaction problems. A 
system of patched and overlaid grids that move relative to each other is used to discretize 
the flow field and the governing equations are integrated using a third-order upwind scheme 
set in an iterative, implicit framework. The code has been used to simulate subsonic flow 
through an axial turbine configuration for which considerable experimental data exists. 
Grid refinement studies have also been conducted as part of the code validation process. 
The current status of the research, along with planned future directions, are also discussed. 
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INTRODUCTION 


• THE PREDICTION OF FLOWS IN TURBOMACHINERY IS A PROBLEM 
OF CONSIDERABLE INTEREST 


• THE FLOWFIELDS ARE INHERENTLY UNSTEADY DUE TO : 

MOTION OF ROTORS RELATIVE TO THE STATORS 

INTERACTION OF DOWNSTREAM AIRFOILS WITH WAKES 
OF UPSTREAM AIRFOILS 

VORTEX SHEDDING AT BLUNT TRAILING EDGES 


• UNSTEADY (OR INTERACTION) EFFECTS PRESENT EVEN FOR AXIAL 
GAPS USED IN CURRENT DESIGNS (50% TO 25% CHORD) 
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INTRODUCTION CONTD. 


• COMPLEX BLADE GEOMETRIES, PERIODIC TRANSITION FROM LAMINAR 
TO TURBULENT FLOW, FLOW SEPARATION, ETC. PRESENT ADDITIONAL 
COMPUTATIONAL DIFFICULTIES 


NEED EXISTS FOR A COMPUTATIONAL TOOL THAT TREATS THE COMPLETE 
TURBINE STAGE AS A SYSTEM AND PROVIDES INFORMATION REGARDING 
THE UNSTEADY EFFECTS 
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POSSIBLE GAINS 


• IMPROVED PERFORMANCE (LOWER LOSSES, ETC.) 


• IMPROVED STRUCTURAL DESIGN (STRUCTURES THAT BETTER WITHSTAND 
THE DYNAMIC STRESSES ASSOCIATED WITH THE UNSTEADY NATURE 
OF THE FLOW) 


• AN UNDERSTANDING OF THE ACOUSTICS (MINIMIZATION OF ENGINE NOISE) 


• PREDICTION OF UNSTEADY HEAT TRANSFER (IMPROVED COOLING TECHNOLOGY) 
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CURRENT STATUS 


• COMPUTER CODES CAPABLE OF PREDICTING UNSTEADY FLOWFIELDS 
IN SINGLE-STAGE TURBOMACHINERY DEVELOPED 

ROTOR-2 2D, NAVIER— STOKES CODE 

ROTOR-3 3D, NAVIER-STOKES CODE 
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CURRENT STATUS CONTD. 


• CODE CAPABILITIES DEMONSTRATED BY SIMULATING AN EXPERIMENTAL 
SINGLE-STAGE AXIAL TURBINE CONFIGURATION 

2D COMPUTATIONS 

3D COMPUTATIONS FULLY ACCOUNTING FOR HUB, CASING AND 
ROTOR TIP CLEARANCE EFFECTS 

GRID— REFINEMENT STUDIES IN 3D 


• STATE-OF-THE-ART COMPUTER GRAPHICS TECHNIQUES USED FOR 
VISUALIZATION OF COMPUTED RESULTS 


"STEREO ANIMATION" TECHNIQUES 


COMPUTATIONAL DETAILS 


HIGH-ORDER UPWIND FINrTE— DIFFERENCE ALGORTHM FOR TIME-ACCURATE 
SOLUTION OF THE THIN-LAYER NAVIER-STOKES EQUATIONS 

ALGORITHM SET IN AN ITERATIVE, FACTORED AND IMPLICIT FRAMEWORK 

SYSTEM OF PATCHED AND OVERLAID GRIDS FOR FLOWFIELD DISCRETIZATION 

ROTOR GRIDS MOVE RELATIVE TO THE STATOR GRIDS 


BALDWIN-LOMAX TYPE TURBULENCE MODEL USED 


COMPUTATIONAL GRID 



ZONE 2 
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SPANWISE VARIATION OF STATOR TIME-AVERAGED PRESSURE 
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PRESSURE AMPLITUDE DISTRIBUTION ON THE STATOR 



X 
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TIME-AVERAGED LIMITING STREAMLINES 
ROTOR SUCTION SURFACE 


FINE GRID 


EXPERIMENTAL VISUALIZATION 
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SUMMARY AND FUTURE DIRECTIONS 


• 2D AND 3D UNSTEADY, THIN-LAYER NAVIER— STOKES CODES DEVELOPED 
TO SIMULATE ROTOR/STATOR INTERACTION IN TURBOMACHINERY 

• CURRENT RESEARCH EFFORTS AND FUTURE DIRECTIONS INCLUDE : 

MULTI-BLADE COMPUTATIONS 
MULTI-STAGE COMPUTATIONS 
BOUNDARY CONDITIONS 
TURBULENCE MODEUNG 
CODE ACCELERATION TECHNIQUES 


GRAPHICS 
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Turbine Stage Aerodynamics 
And Heat Transfer Prediction 


L. W. Griffin and H. V. McConnaughey 
NASA/MSFC 


A numerical study of the aerodynamic and thermal environment 
associated with axial turbine stages is presented. 
Computations were performed using a modification of the 
unsteady NASA Ames viscous code, R0T0R1, and an improved 
version of the NASA Lewis steady inviscid cascade system 
MERIDL-TSONIC coupled with boundary layer codes BLAYER and 
STAN5 . Two different turbine stages were analyzed: the 
first stage of the United Technologies Research Center Large 
Scale Rotating Rig (LSRR) and the first stage of the Space 
Shuttle Main Engine (SSME) high pressure fuel turbopump 
turbine. The time-averaged airfoil midspan pressure and heat 
transfer profiles were predicted for numerous thermal 
boundary conditions including adiabatic wall, prescribed 
surface temperature, and prescribed heat flux. Computed 
solutions are compared with each other and with experimental 
data in the case of the LSRR calculations. Modified R0T0R1 
predictions of unsteady pressure envelopes and instantaneous 
contour plots are also presented for the SSME geometry. 
Relative merits of the two computational approaches are 
discussed . 
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OVERVIEW 


• OBJECTIVES 

• MOTIVATION 

• APPROACH 

• CODES 

• APPLICATIONS 

• RESULTS 

• AERODYNAMIC 

• HEAT TRANSFER 

• COMPUTER RESOURCES REQUIREMENTS 

• IMPACT 

• SUMMARY AND CONCLUSIONS 


5 - 1556 - 9-105 


219 


OBJECTIVES 


• DEFINITION OF STEADY AND UNSTEADY AEROTHERMAL 
ENVIRONMENTS FOR THE SPACE SHUTTLE MAIN ENGINE (SSME) 
HIGH PRESSURE FUEL TURBOPUMP (HPFTP) TURBINE 

• CALIBRATION OF TURBINE ANALYSIS CODES VIA COMPARISON 
OF PREDICTED AND EXPERIMENTALLY MEASURED AIRFOIL 
LOADS IN UNITED TECHNOLOGIES RESEARCH CENTER LARGE 
SCALE ROTATING RIG (LSRR) 


5 - 1557 - 9-105 
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MOTIVATION 
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APPROACH 


CODES: 

• MTSBL- - 
STAN5 

• ROTOR1 - 

APPLICATIONS: 

• LSRR 

• SSMEHPFTP 


STEADY INVISCID QUASI-3D SYSTEM MERIDL- 
TSONIC COUPLED WITH BOUNDARY LAYER CODES 
BLAYER AND STAN5. DEVELOPED AT NASA/LEWIS. 

UNSTEADY 2D SINGLE-STAGE CODE DEVELOPED AT 
NASA/AMES. ASSUMES EQUAL BLADE COUNT. 
MODIFIED TO INCORPORATE HEAT TRANSFER 
PREDICTION CAPABILITY 


- CODE CALIBRATION 

- ENVIRONMENT DEFINITION 


5 - 1558 - 9-105 
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PRESSURE DISTRIBUTIONS FOR THE LSRR 




Stator Rotor 
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PRESSURE DISTRIBUTIONS FOR 
THE SSME HPFTP 




Stator 


Rotor 
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UNSTEADY PRESSURE DISTRIBUTIONS ON 
THE SSME HPFTP ROTOR 



x/c 


5 - 1 592 - 9-1 05 


225 


ENTROPY CONTOURS 
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COMPARISON OF PREDICTED AND EXPERIMENTAL 
HEAT TRANSFER FOR THE LSRR 




Stator Rotor 


5 - 1564 - 9-105 
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COMPUTER RESOURCES REQUIREMENTS 


MTSBL - 
ROTOR1 


STAN5 - 6 CPU MINUTES, 0.9 X 10 6 WORDS MEMORY 

-6 CPU HOURS, 2.1 X10 6 WORDS MEMORY + 

4.1 X10 6 WORDS SSD 


5 - 1559 - 9-105 
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IMPACT 


• MSFC BENCHMARKING OF QUASI-3D TURBINE ANALYSIS SYSTEM 
AND UNSTEADY STAGE CODE; ASSESSMENT OF RELATIVE 
MERITS OF THE TWO APPROACHES. 

• STEADY 3D FIRST AND SECOND STAGE BLADE LOADS WERE 
USED IN STRESS ANALYSIS OF HPFTP TURBINE BLADES IN 
SUPPORT OF RETURN TO FLIGHT : INTEGRATED LOADS PREDICT 
CORRECT TORQUE WITHIN 1%. 

• UNSTEADY MIDSPAN PRESSURES WERE USED IN STRUCTURAL 
DYNAMICS ANALYSIS OF HPFTP FIRST STAGE BLADES. 
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SUMMARY AND CONCLUSION 


• SUPPORTED MSFC SSME PROGRAMS THROUGH 
CHARACTERIZATION OF THE STEADY AND UNSTEADY 
AEROTHERMAL ENVIRONMENTS OF THE HP FTP. 

• BOTH APPROACHES YIELDED COMPARABLE STEADY OR TIME- 
AVERAGED AERODYNAMIC AND THERMAL RESULTS. 

• AERODYNAMIC RESULTS FROM BOTH APPROACHES COMPARE 
WELL WITH DATA; THERMAL RESULTS EXHIBIT FAIR AGREEMENT 
WITH DATA. 
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AUTOMATED DESIGN OF CONTROLLED DIFFUSION BLADES 


Jose M. Sanz 

NASA Lewis Research Center 


ABSTRACT 


A numerical automation procedure has been developed to be 
used in conjuction with an inverse hodograph method for the 
design of controlled diffusion blades. 

With this procedure a cascade of airfoils with a prescribed 
solidity , inlet Mach No., inlet air flow angle and air 
flow turning can be produced automatically. 

The trailing edge thickness of the airfoil, an important 
£ quantity in inverse methods, is also prescribed. 

The automation procedure consists of a multi-dimensional 
Newton iteration in which the objective design conditions are 
achieved by acting on the hodograph input parameters of the 
underlying inverse code. 

The method, although more general in scope, is applied 
in this paper to the design of axial flow turbomachinery 
blade sections, both compressors and turbines. A collaborative 
effort with U.S. Engine Companies to identify designs of 
interest to the industry will be described. 
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CURRENT APPROACH 


• GUESS GEOMETRY 

- SINGLE CIRCULAR ARC 

- MULTIPLE CIRCULAR ARC 

- POLYNOMIAL SHAPES 

- NACA AIRFOILS SERIES 


• ANALYZE SHAPE WITH FLOW SOLVER 


- 2-D, QUASI 3-D OR 3-D CODES 

- POTENTIAL, EULER OR NAVIER-STOKES 
SOLVERS 


BUILDS IN EMPIRICISM 


AUTOMATION PROCEDURE 


ITERATION OF INVERSE HODOGRAPH METHOD 


PRODUCES BLADE WITH PRESCRIBED 

- SOLIDITY 

- INLET MACH NUMBER 

- INLET AIR ANGLE 

- AIR FLOW TURNING 

- TRAILING EDGE THICKNESS 
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BENEFITS 


• REPLACE TRIAL AND ERROR PROCESS BY 
EXACT 2-D SOLUTION 

• GREATLY REDUCES MANPOWER AND 
TURN AROUND TIME 

9 INITIALLY OPTIMIZED BLADE SECTIONS 

• INNOVATIVE DESIGNS 


MAJOR IMPACT FOR DESIGNS PROBLEMS WITH 
NO DATA BASE 


TO IMPLEMENT AUTOMATION 


EQUATION U (R, M Q , 9, q (s)) 

IS REPLACED BY 

y = F (X) 


WHERE 


X = (R, M 0 , 0, Q m , Q te , S te ) 


y = y (o, M v (3 r Ap, dn te * ds te ) 


CD-88-33796 
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NEWTON ITERATION TO SOLVE VECTOR EQUATION 


F(x) - y 0 = o, 

y 0 = (o, M r a ( 3< dn te , ds te ) 0 

IS OBJECTIVE FUNCTION. 

ITERATION 

x n+1 x n " ^ 

JACOBIAN: J = / dF j 

\ dx j 
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APPLICATIONS 


• ARC 40X80X120 WIND TUNNEL TURNING VANES 


• SUBSONIC AND TRANSONIC (SHOCK-FREE) 

TURBOMACHINERY AND PROPELLER SECTIONS 


• INLET GUIDE VANES FOR LARGE RANGE OF 
INLET AIR FLOW ANGLE - EXPERIMENTAL 
VERIFICATION 3/88 


ARC 40X80X120 WIND TUNNEL TURNING VANES 
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AUTOMATED DESIGN OF TURBINE BLADE SECTIONS 
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AUTOMATED DESIGN OF TURBINE BLADE SECTIONS 
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Numerical Analysis of 
Cascade 


Flow Through Oscillating 
Sections 


by Dennis L. Huff 
NASA Lewis Research Center 


The design of turbomachinery blades requires the prevention 
of flutter for all operating conditions. However, flow field 
predictions used for aeroelastic analysis are not well understood 
for all flow regimes. The present research focuses on numerical 
solutions of the Euler and Navier-Stokes equations using an ADI 
procedure to model two-dimensional, transonic flow through 
oscillating cascades. The model prescribes harmonic pitching 
motions for the blade sections for both zero and non-zero 
inter-blade phase angles. The code introduces the use of a 
deforming grid technique for convenient specification of the 
periodic boundary conditions. Approximate nonreflecting boundary 
conditions have been coded for the inlet and exit boundary 
conditions. Sample unsteady solutions have been performed for an 
oscillating cascade and compared to experimental data. Also, 
test cases were run for a flat plate cascade to compare with an 
unsteady, small-perturbation, subsonic analysis. 

The predictions for oscillating cascades with non-zero 
inter-blade phase angles are in good agreement with experimental 
data and small-perturbation theory. The zero degree inter-blade 
phase angle cases, which were near a resonant condition, differ 
from the experiment and theory. Studies on reflecting versus 
non-reflecting inlet and exit boundary conditions show that the 
treatment of the boundary can have a significant effect on the 
first harmonic, unsteady pressure distributions for certain flow 
conditions. This code is expected to be used as a tool for 
reviewing simpler models that do not include the full non-linear 
aerodynamics or as a final check for designs against flutter in 
turbomachinery . 
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OBJECTIVE: 

• DEVELOP AN UNSTEADY, VISCOUS, TRANSONIC FLOW ANALYSIS METHOD 
FOR CASCADED AIRFOILS WITH PITCHING MOTIONS 

TECHNIQUE: 

• A COMPRESSIBLE, FULL NAVIER-STOKES, FINITE DIFFERENCE CODE IS 
DEVELOPED TO MODEL OSCILLATING CASCADES FOR BOTH ZERO AND NON- 
ZERO INTER-BLADE PHASE ANGLES 

• A UNIQUE DEFORMING GRID TECHNIQUE IS INTRODUCED TO CAPTURE 
BLADE MOTIONS 

• APPROXIMATE NON-REFLECTING BOUNDARY CONDITIONS ARE USED TO 
MINIMIZE WAVE REFLECTIONS FROM THE OUTER BOUNDARIES 


CD-88-38197 
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DEFORMING GRID TECHNIQUE, SIMPLIFIED GRID WITH 
EXAGGERATED MOTION 



\ 
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BOUNDARY CONDITIONS 

INLET: 

• SPECIFY DENSITY, VELOCITY, AND FLOW ANGLE 

• SOLVE THE CHARACTERISTIC FORM OF THE GOVERNING EQUATIONS TO 
DETERMINE THE ENERGY 


EXIT: 

• FOR VISCOUS FLOWS, SPECIFY STATIC PRESSURE AND EXTRAPOLATE 
DENSITY AND VELOCITY FROM THE INTERIOR 

• FOR INVISCID FLOWS, SPECIFY STATIC PRESSURE AND SOLVE THE 
CHARACTERISTIC FORM OF THE GOVERNING EQUATIONS TO DETERMINE 
DENSITY AND VELOCITY 

PERIODIC BOUNDARIES FROM BLADE-TO-BLADE 

SOLID WALL BOUNDARIES ON AIRFOIL SURFACE 

AVERAGE FLOW PROPERTIES ACROSS THE WAKE 
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RESULTS 


1. NACA 65-SERIES CASCADE, M = 0.17, TWO INTER-BLADE PHASE ANGLES, 
COMPARISONS WITH EXPERIMENT 

2. FLAT PLATE CASCADE, M = 0.65 AND 0.80, THREE INTER-BLADE PHASE 
ANGLES, COMPARISONS WITH SMALL-PERTURBATION THEORY 

3. NASA LEWIS CASCADE, M = 0.65 AND 0.80, THREE INTER-BLADE PHASE 
ANGLES, BOUNDARY CONDITION STUDY, NUMERICAL TIME ACCURACY STUDY, 
COMPARISONS WITH EXPERIMENT. 
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MEAN FLOW PRESSURE DISTRIBUTION FOR 
NACA 65-SERIES CASCADE 

M, = 0.17, y = 55°, g/c = 0.75, r' = 0.06 
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FIRST HARMONIC PRESSURE DISTRIBUTION FOR 
NACA 65-SERIES OSCILLATING CASCADE 

M-) = 0.17, 0i = 64° ±2°, k = 0.123 
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BOUNDARY CONDITION AND THEORY COMPARISON FOR 
FLAT PLATE CASCADE 

M 1 = 0.65, 0 -\ = 530 ±0.100, k = 0.221 



x'/c 
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BOUNDARY CONDITION AND THEORY COMPARISON FOR 

FLAT PLATE CASCADE 

M 1 = 0.65, /?! = 53° ±0.10°, k = 0.221 
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MEAN FLOW PRESSURE DISTRIBUTION FOR BICONVEX 

AIRFOIL CASCADE 

/ 3 1 = 60°, 7 = 53°, g/c = 0.767, r' = 0.07 



x'/c 
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FIRST HARMONIC PRESSURE DISTRIBUTION FOR 
BICONVEX AIRFOIL CASCADE 

M, = 0.80, 0, = 60°±1.2®, k = 0.185. o = -90° 






CD-88-38219 


256 


FIRST HARMONIC PRESSURE DISTRIBUTION FOR 
BICONVEX AIRFOIL CASCADE 

M-, = 0.80, 0i = 60° ±1.2°, k = 0.185, a = 0° 
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ANALYSIS, UPPER SURFACE 
ANALYSIS, LOWER SURFACE 
EXPERIMENT, UPPER SURFACE 
EXPERIMENT, LOWER SURFACE 
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CONCLUSIONS 


• A COMPRESSIBLE EULER OR NAVIER-STOKES, FINITE-DIFFERENCE CODE HAS BEEN 
DEVELOPED FOR THE ANALYSIS OF OSCILLATING CASCADES 

• A DEFORMING GRID TECHNIQUE IS USED TO CAPTURE THE MOTION OF THE AIRFOILS 
FOR BOTH ZERO AND NON-ZERO INTER-BLADE PHASE ANGLES 

• APPROXIMATE TWO-DIMENSIONAL, UNSTEADY CHARACTERISTIC BOUNDARY CONDITIONS 
ARE USED AT THE INLET AND EXIT TO MINIMIZE WAVE REFLECTIONS 

• IN GENERAL, PREDICTIONS WITH NON-ZERO- INTER-BLADE PHASE ANGLES ARE IN GOOD 
AGREEMENT WITH THE EXPERIMENTAL DATA AND SMALL-PERTURBATION THEORY 

• PREDICTIONS FOR ZERO DEGREE INTER-BLADE PHASE ANGLE CASES, WHICH WERE NEAR 
AN ACOUSTIC RESONANT CONDITION, DIFFER FROM THE EXPERIMENT AND THEORY 

• THE TYPE OF BOUNDARY CONDITION USED AT THE INLET AND EXIT CAN HAVE A SIG- 
NIFICANT EFFECT ON THE FIRST HARMONIC, UNSTEADY PRESSURE DISTRIBUTIONS 

• FIRST-ORDER AND SECOND-ORDER TEMPORAL ACCURACY RESULTS DO NOT SHOW SIGNIF- 
ICANT DIFFERENCES IN THE UNSTEADY PRESSURE DISTRIBUTIONS 
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ANALYSIS OF THREE-DIMENSIONAL VISCOUS FLOW 
IN A SUPERSONIC THROUGHFLOW FAN 

by 

Dr. Rodrick V. Chima 
NASA Lewis Research Center 
Cleveland, OH 44135 


Abstract 

A three-dimensional Navier-Stokes code has been developed for analysis of turboma- 
chinery blade rows and other internal flows. The Navier-Stokes equations are written in 
a Cartesian coordinate system rotating about the x-axis, and then mapped to a general 
body-fitted coordinate system. Streamwise viscous terms are neglected using the thin- 
layer assumption, and turbulence effects are modelled using the Baldwin-Lomax turbu- 
lence model. The equations are discretized using finite differences on stacked C-type grids 
and are solved using a multistage Runge-Kutta algorithm with a spatially-varying time 
step and implicit residual smoothing. 

Calculations were made of the flow around a supersonic throughflow fan blade. The 
fan was designed at NASA Lewis Research Center as a key component in a supersonic 
cruise engine. It was designed to produce a total pressure ratio of 2.7 at an axial Mach 
number of 2.0. The midspan section of the blade is being tested in a supersonic linear 
cascade at Virginia Polytechnic Institute and will be tested in a rotating rig at Lewis in 
the near future. Comparisons between earlier quasi-3-D calculations and the VPI data 
show excellent agreement between shock locations and wake traverses. 

The 3-D calculations were done on a 129 * 29 * 33 grid and took 50 minutes of cpu 
time on a Cray X-MP. Comparisons with the quasi-3-D results show minor differences in 
loading due to 3-D effects. Particle traces show nearly 2-D flows near the pressure surface, 
but large secondary flows within the suction surface boundary layer. The horseshoe vortex 
ahead of the leading edge is clearly seen. 

References 

1. Chima, R. V., and Yokota, J. W. “Numerical Analysis of Three- Dimensional Viscous 

Internal Flows,” NASA TM-100878, July, 1988. 

2. Schmidt, J. F., Moore, R. D., and Wood, J. R. “Supersonic Throughflow Fan Design,” 

NASA TM-88908, AIAA-87-1746, June, 1987. 
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RVC3D (ROTOR VISCOUS CODE 3-D) 

BY R. V. CHIMA 
DESCRIPTION 

• EULER OR NAVIER-STOKES ANALYSIS FOR STEADY 
3-D FLOWS IN TURBOMACI11NERY BLADE PASSAGES 

FEATURES 

• STACKED C-TYPE GRIDS FOR AXIAL OR CENTRIFUGAL MACHINES 

• CARTESIAN FORMULATION ROTATING ABOUT X-AXIS 
RECTANGULAR OR ANNULAR GEOMETRIES 

• SOLVES NAVIER-STOKES EQUATIONS IN FINITE-DIFFERENCE FORM 
THIN-LAYER FORMULATION NEGLECTS STREAMWISE VISCOUS TERMS 
RETAINS IIUB-TO-TIP k BLADE-TO-BLADE VISCOUS TERMS 
BALDWIN-LOMAX TURBULENCE MODEL 

• EXPLICIT 4-STAGE RUNGE-KUTTA TIME-MAllCHING SCHEME 
VARIABLE At itj k IMPLICIT RESIDUAL SMOOTHING 
HIGHLY VECTORIZED FOR CRAY X-MP 


RESULTS 


. SUPERSONIC THROUGHFLOW FAN 


GOVERNING EQUATIONS 


d t q + J[dtE + d n F + d s G- Re~ 1 (d r) Fv +d i G v )] = H 
WHERE: 

q = [p, pu , pv , pw, e] T 

H— [0, 0, —Upw, O pv, 0] T 

E = J _1 [pU 1 , puU' + £xP, pvU' + £ y p, pu)U' + izP, eU' +pU] T 
F = J~ l [ pV', puV' + q x p, pvV' + T] y p, pwV 1 + TjzP, eV' + pV f 
G = J- 1 [ pW', puW' + f s p, pvW' + f„p, pwW' + &p, eW' + P W f 

RELATIVE VELOCITIES: 

Xl! — u 
v = v — Clz 
w' = w + 0 y 

RELATIVE CONTRAVARIANT VELOCITIES: 

U' = i x u + i y v' + i z w' 

V' = T] x u + r) y v' + rizW 1 
W' = $X U + SyV 1 + £gli/ 

ENERGY AND STATIC PRESSURE: 
e = p [C V T + (u 2 + v 2 + it> 2 )/2] 
p = (7 - 1) [e - p(tt 2 + v 2 + w 2 )/2] 
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MULTISTAGE RUNGE-KUTTA ALGORITHM 


GOVERNING EQUATIONS 

dtq = —J [72/ — ( Rv + 77)] 

72/ = INVISCID RESIDUAL 
Rv = VISCOUS RESIDUAL 
D = ARTIFICIAL DISSIPATION TERM 
MULTISTAGE SCHEME 


<70 = Qn 


qi — qo — ot\J At [72/ q 0 - (R v + D) g 0 ] 


qk-qo- a k J At [72/ q k — i - (72^ + D) q 0 ] 


qn+i — qk 


Rv & D EVALUATED AT FIRST STAGE ONLY 
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ARTIFICIAL DISSIPATION 


NONCONSERVATIVE VERSION OF JAMESON FORMULATION 
Dq = + Dtj + D$) q 

^-DIRECTION OPERATOR 
Dzq = Ct [V 2 q^ - V 4 q^ a ) 

WHERE: 

q = I(_L + _L L!(J_ + ±) 

{ AtJ J\A Srl AsJ 

V 2 = V 2 max [v i+ 1 , Vi , //,•_ i ) 

V 4 = max (0, fj. 4 — V 2 ) 

AND 

. _ lA'+i.y ~ 2-Pt,; + Pj—\,j | 
l,; l-Pt+i,; + 2F t -,y + -Pt-i.yl 


/*2 = 0(1) 
/*4 = 0(i) 
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IMPLICIT RESIDUAL SMOOTHING 


USE A TIME STEP GREATER THAN THE STABILITY LIMIT 
MAINTAIN STABILITY BY SMOOTHING THE RESIDUAL IMPLICITLY 


{1 - e^6 u )(l - - e^)R = R 

UNCONDITIONALLY STABLE IF 



WHERE 

A* IS COURANT LIMIT OF THE UNSMOOTHED SCHEME 
A IS THE LARGER OPERATING COURANT NUMBER 
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SUPERSONIC FAN 

(NASA LeRC ASSESSMENT STUDY) 


SUPERSONIC 

DIFFUSER-.. 


SUBSONIC 

DIFFUSER 


,-3 STAGE CONVENTIONAL FAN 





r -SINGLE STAGE SUPERSONIC FAN 
\ r- DUCT NOZZLE 

4 - 


CONVENTIONAL TURBOFAN 


SUPERSONIC , 
DIFFUSER 



SUPERSONIC FAN ENGINE 


SUPERSONIC FAN ENGINE FEATURES 

• SHORT, ALL SUPERSONIC INLET 

• SINGLE STAGE SUPERSONIC FAN 

• BPR DECREASES WITH Mq 


IMPLICATIONS 

• LOWER WEIGHT, LOWER INLET DRAG 

• LOWER WEIGHT AND COST, RUGGED BLADING 

• HIGHER CRUISE THRUST 


WEIGHT 15 TO 20% LOWER CRUISE SFC-10 TO 20% LOWER 




GEOME-T^ 

thr oug«flcw FM< 
SUPERSONIC THRO 




pressure 


SUPERSONIC THROUGHFLOW FAN. M,.l » 2-5 


to 

ON 

-4 
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DENSITY 


SUPERSONIC THROUGH FLOW FAN, Mr e 1 = 2.5 
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PARTICLE TRACES 


SUPERSONIC THROUGHFLOW FAN, Mr el = 2.5 
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PARTICLE TRACES 
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SUMMARY 


PHYSICS 

. 3-D NAVIER-STOKES ANALYSIS FOR STEADY INTERNAL FLOWS 

• CARTESIAN FORMULATION ROTATING ABOUT X-AXIS 

• THIN LAYER IN STREAMWISE DIRECTION, FULL N-S IN OTHERS 

• BALDWIN-LOMAX TURBULENCE MODEL 

NUMERICS 

• FINITE-DIFFERENCE FORM ON GENERAL BODY FITTED GRID 

• EXPLICIT MULTISTAGE RUNGE-KUTTA SCHEME 

• VARIABLE A t itj L IMPLICIT RESIDUAL SMOOTHING 

RESULTS 

. SUPERSONIC THROUGHFLOW FAN 

• HORSESHOE VORTEX AHEAD OF CYLINDER 

• ANNULAR TURBINE CASCADE 

FUTURE 

. NEW FAN DESIGN WITH CONVERGING HUB, HIGH TURNING 

• FINER GRIDS ON NAS 


• TIP CLEARANCE 

• MULTIGRID 


SESSION X 


STOVL 


Chairman: 

Meng-Sing Liou 

Chief, Computational Fluid Dynamics Branch 
NASA Lewis Research Center 
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SIMULATION OF POWERED-LIFT FLOWS 


William R. Van Dalsem 
Kalpana Chaw la 
Karlin R. Roth 
Merritt H. Smith 
Kuditipudi V. Rao 

NASA Ames Research Center, Moffett Field, CA 
Thomas C. Blum 

BOEING Advanced Systems, Seattle, WA 


The primary objective of this presentation is to expose government, industry, and 
academic scientists to work underway at NASA- Ames towards the application of CFD to 
the powered lift area. One goal of our research is to produce the technologies which will be 
required in the application of numerical techniques to, for example, the Supersonic STOVL 
program. 

In the presentation, we will summarize our progress to date on the following specific 
projects: 

Jet in ground effect with crossflow 
Jet in a crossflow 

Delta planform with multiple jets in ground effect 
Integration of CFD with thermal and acoustic analyses 
Improved flow visualization techniques for unsteady flows 
YAV-8B Harrier simulation program 
E-7 simulation program 

Additional work is underway at NAS A- Ames in the development of turbulence models and 
solution adaptive grid techniques suitable for the powered lift area, and the simulation of 
USB configurations. However, this work is not included here due to space constraints. 
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POWERED-LIFT CFD PROGRAM 

OBJECTIVE 

Develop and validate CFD technologies for the V /STOL 
field, with particular emphasis on the requirements 
of the supersonic STOVL program. 

TECHNICAL APPROACH 

Time-accurate Navier- Stokes solutions on overlapped 
adapting grids, coupled to thermal, propulsive, con- 
trols, and acoustic analyses. 

STRATEGY 

• Starting from the RFA CFD technology base, develop in-depth expertise in the most 
critical components of the flow about V/STOL aircraft. 

• As required, improve the computational techniques to allow the accurate simulation 
of these flows. 

• Simultaneously, begin complete V/STOL aircraft simulation efforts (e.g., VSRA 
and E-7), with above “component” efforts contributing required technologies and 
expertise. 



IN A TIMELY MANNER, DEVELOP BOTH IN-DEPTH EXPERTISE AND THE 
ABILITY TO TREAT COMPLETE V/STOL CONFIGURATIONS. 


ORIGINAL PAGE IS 

OF POOR QUALITY 


LANDING SURFACE I VEHICLE 


POWERED-LIFT CHALLENGES & CFD REQUIREMENTS 



FUNDAMENTAL FLUID MECHANICS 


JETS EXITING 

IMPACT OF 

MULTIPLE JETS/ 



FROM A 

BODY DYNAMICS 

FOUNTAIN 

EJECTORS 


SURFACE 


EFFECTS 




POWERED-LIFT DESIGN PROBLEMS 

SUCK DOWN 
FOUNTAIN EFFECTS 

LIDS AND LOWER SURFACE CONTOURS 
GROUND VORTEX 

INDUCED LIFT AND PITCHING MOMENTS 
LATERAL AND DIRECTIONAL CHARACTERISTICS 
HOT GAS INGESTION 

LANDING GEAR AND LANDING ENVIRONMENT 
LOWER SURFACE HEATING 
ACOUSTIC LOADS 
ETC ... 


< COMPLEX 
GEOMETRIES 

SOLUTION 

ADAPTIVE 

MOVING GRID 
TECHNOLOGY 

UNSTEADY 

GRAPHICS 

GRID GEI 

NERATION 

DYNAMICS 


CFD TECHNOLOGY REQUIREMENTS 
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JET IN GROUND EFFECT WITH CROSSFLOW 


A fundamental component of STOVL take-off/landing flow physics 

• To date, CFD has proved capable of resolving the salient features, and 
adding to previous understanding: 

O Computed correct ground vortex formation and extent over a 
range of jet to freestream velocity ratios and jet heights. 

O Allowed a systematic study of the impact of a variety of flow conditions, 
including jet shape, moving ground board, etc... 

° Simulated ring vortex and shock disc motion, which may be important 
sources of unexplained intense noise levels 

• In the future, will focus on improving our understanding of jet unsteadiness and 
the existing discrepancies between experimental and full-scale studies. 



IMPACT OF TURBULENCE MODELING 

C ALONG CENTERLINE OF JET 
P 

V e - 0.223. h/D - 3 
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JET IN GROUND EFFECT EMULATION 

(Ve=0.10, h/D=3, IrbHnt i«t) 
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UNSTEADY FLOW ANALYSIS 

BACKGROUND The unsteady flow about powered-lift vehicles induces significant 

unsteady loads, hence must be accurately predicted. 

Require significant software improvements to deal with unsteady 
flows on a routine basis. 

APPROACH Develop techniques which allow the rapid yet indepth analysis and 

validation of unsteady flow simulations. 

PL0T4D 

An “upgraded” version of PL0T3D program which allows 
the interactive animation and analysis of unsteady data sets. 



STREAKER 

Produces time-accurate numerical streaklines (smoke) suitable 
for direct comparison with experimental unsteady flow visualizations. 


NUMERICAL STREAKLINES 




' EXPERIMENTAL SMOKE FLOW VISUALIZATION 


NUMtKILML sincMrviLiNLj 
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THERMAL LOADS 


PROVEN RISK OF VEHICLE AND LANDING SURFACE DAMAGE (e.g., LPH TRIPOLI DECK BUCKLING) 
DURING TAKE-OFF/LANDING, AND VEHICLE DAMAGE DURING TRANSITION 

• AT PRESENT, EXPENSIVE FULL-SCALE JET ENGINE/MATERIALS TEST (FLUK) ARE ONLY 
METHOD OF TESTING UNSTEADY THERMAL RESPONSE OF MATERIALS TO ENGINE EXHAUST 

• COUPLED UNSTEADY JET FLOW CFD/THERMAL ANALYSIS OF SURFACES MADE OF A VARIETY 
OF MATERIALS (FROM ALUMINUM AND COMPOSITES, TO CONCRETE) WOULD ALLOW DETAILED 
ANALYSIS OF V/STOL THERMAL DAMAGE 


NAVIER STOKES EOS. 

°t + E x + F y + G z = Re_1 < Ev x + Fv y + Gv z > 




ORIGINAL PAGE IS 
OF POOR QUALITY 




H»T JET THERMAL 

(Ve=0.10, h/D=3, tlrb 


NAVIER-STOKES EQS 
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JET IN A CROSSFLOW 


• Isolates physics of aero/propulsive interaction region during transition. 


• Navier-Stokes solution captures all critical flow features. 


• Quantitative correlation with existing database demonstrated. 


• Current efforts focus on computing realistic STOVL geometries. 


FLOW PHYSICS 



JET CENTERLINE 

Mj/Moo = 4.0, Moo = 0.19 



K. Roth, FFF 


284 ORIGINAL PAGE IS 

OF POOR QUALITY 





285 


GROUND EFFECTS 


C, AND C m OF V/STOL VEHICLES IN GROUND-EFFECT ARE STRONGLY INFLUENCED BY 
HEIGHT ABOVE GROUND AND RATE OF ASCENT/DESCENT 

• DIFFICULT TO STUDY EXPERIMENTALLY, AT PRESENT GOVERNING FLOW PHYSICS IS 
NOT UNDERSTOOD 

• TO DATE, CFD HAS SHED INSIGHT IN TO FLOW PHYSICS OF NEGATIVE C, AT LOW 
HEIGHTS, AND WORK IS IN PROGRESS TO STUDY DYNAMIC CASE 

0 (a, h, V e ) = EFFECTIVE ANGLE OF ATTACK 
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E-7 COMPUTATION THROUGH FULL-SCALE EXPERIMENT PROGRAM 
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SUMMARY 


• To date, have studied the numerical simulation of the following powered-lift flows: 

O Jet in ground effect with crossflow (including thermal loads) 

O Jet in crossflow 

O Delta with multiple jets in ground effect with crossflow 

Comparison with data indicates that these simulations predict the fundamental flow 
phenomena and yield quantitative results for many of the observed trends. 

• This experience has motivated us towards additional R&D in the following areas: 

o Adaptive gridding 
O Improved turbulence modeling 
O Unsteady flow simulation and analysis 

• Presently also working towards simulation of the following "complete" configurations: 


O STOVL CFD Validation model (In cooperation with FFF) 

O YAV-8B Harrier 

( 

it 

FAP to obtain data) 

O E-7 

( 

ii 

Boeing) 

O USB configuration 

( 

ii 

FAP) 

O Tiltrotor configuration 

( 

ii 

RFA Rotorcraft group) 
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A Numerical Study of the Hot Gas Environment 
Around a STOVL Aircraft in Ground Proximity 

hy 

Thomas J . VanOverbeke 
and 

James D. Holdeman 
of 

NASA Lewis Research Center 


ABSTRACT 

The development of Short Take-off Vertical Landing (STOVL) aircraft 
has been based on empiricism. In this study, a 3-D flow code was used to 
calculate the hot gas environment around a STOVL aircraft in ground 
proximity. Preliminary calculations are reported to identify key 
features of the flowfield, and to demonstrate the capability of a CFD 
code to calculate the temperature of the gases ingested at the engine 
inlet for typical flow and geometric conditions. 
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CALCULATIONS WERE DONE WITH A 3-D TEACH-TYPE CODE 
TO SOLVE THE TIME-AVERAGED NAVIER-STOKES EQUATIONS 

EQUATIONS ARE SOLVED SEQUENTIALLY 

A PRESSURE CORRECTION EQUATION IS SOLVED USING 
THE SIMPLE ALGORITM OF PATANKAR 

HYBRID DIFFERENCING WAS USED TO SAVE MEMORY 

A TWO EQUATION TURBULENCE MODEL WAS USED 

CODE WAS USED FOR COMBUSTORS AND DID NOT INCLUDE 
A DENSITY CORRECTION TERM FOR RAPIDLY VARYING 
PRESSURE 

EXIT BOUNDARY CONDITIONS WERE CHANGED FOR 
MODELING THE HOT GAS INGESTION PROBLEM 
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NASA 

C-88-06794 


SCHEMATIC OF COMPUTATIONAL DOMAIN 
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NEAR-GROUND PLANE VELOCITY AND TEMPERATURE DISTRIBUTIONS 

(H/Dj = 4) 



NEAR-GROUND PLANE VELOCITY AND TEMPERATURE DISTRIBUTIONS 

(H/Dj=4) 
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CALCULATED VELOCITY AND TEMPERATURE DISTRIBUTIONS 
IN VERTICAL X-Y PLANE THROUGH ENGINE INLET 

(H/Dj=4; lloo/Vj=0.03) 



TAIN 


A) VELOCITY 
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B) TEMPERATURE 
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TEMPERATURE DISTRIBUTIONS IN VERTICAL PLANES FROM 
AIRCRAFT CENTERPLANE TO OUTBOARD OF FUSELAGE 

(H/Dj = 4) 
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Uoo/Vj = 0.03 


1100 °F 


0.09 
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TEMPERATURE DISTRIBUTIONS IN VERTICAL PLANES 
FROM AIRCRAFT CENTERPLANE TO OUTBOARD OF FUSELAGE 

(U°o/Vj=0.09) 
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TEMPERATURE DISTRIBUTION IN HORIZONTAL PLANES 
EFFECT OF CHANGING HEIGHT 


1100 °F 


0 °F 



HEADWIND 



Uoo/Vj = 0.09 



H/Dj = 4 





H/Dj = 2 
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TEMPERATURE DISTRIBUTIONS IN VERTICAL 
(Y-Z) PLANES THROUGH ENGINE INLET 


Uoo/Vj=0.03 0.09 



C ni 
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Jet, Ambient, and Inlet 
Temperatures for Cases Calculated 3 


UjVj 

H/Dj 

u m , 

kn 

T 

1 avg 

T 

1 max 

T • 

1 nun 

0.03 

4 

17.8 

187.7 

388.2 

109.8 

.03 

2 

17.8 

173.5 

429.6 

72.0 

.09 

4 

53.6 

145.1 

469.9 

61.6 

.09 

2 

53.6 

97.6 

375.3 

62.3 


a At nozzle pressure ratio of 1.21; 7} = 1000; = 70 

u> (all lemperatures arc degrees F). 

o 

v© 


Average Inlet Temperatures 3 


Distance from exhaust 

Ratio of forward speed of 

lift jets to 

aircraft (or strength of 

ground plane, 

headwind) to exhaust 

H/Dj 

jet velocity, 


u m 



0.03 

0.09 

4 

0.13 

0.08 

2 

.11 

.03 


a (A,, - TJHJj - TJ 


Inlet Temperature Distortion 3 


Distance from exhaust 

Ratio of forward speed of 

lift jets to 

aircraft (or strength of 

ground plane, 

headwind) to exhaust 

H/Dj 

jet velocity, 



/Vj 


0.03 

0.09 

4 

0.30 

0.44 

2 

.38 

.34 


a (T™ - T m ,nV(Tj - Tj 
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SUMMARY 


INTERNAL FLOW CODE CAN BE USED IN PREDICTING STOVL FLOWFIELDS 
SIGNIFICANT INGESTION WAS PREDICTED IN ALL CASES CALCULATED 
ALL CASES PREDICTED NEAR FIELD INGESTION 


WEAKER HEADWINDS ALSO PREDICTED INGESTION BY GROUND VORTEX FLOW 
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CFD ANALYSIS FOR HIGH SPEED INLETS 


Tom Benson 


NASA Lewis Research Center 


ABSTRACT 


The increased national interest in high speed flight has increased research for 
high speed propulsion components. The highly three-dimensional flows present in 
supersonic/hypersonic inlets are currently being studied at NASA Lewis both 
experimentally and computationally using a family of steady PNS and NS solvers 
and unsteady NS solvers. This paper presents some of the results of these efforts 
with an emphasis on the comparison of the computational and experimental 
results. 

The flow in high speed inlets typically involves the interaction of compression 
shock waves and boundary layers on the internal surfaces. The fundamentals of 
these interactions have been studied experimentally for many years, while more 
recently, computations have been used to study these complex three dimensional 
flow fields. Attempts to control the flow through boundary layer bleed are being 
investigated computationally prior to wind tunnel experiments. The ultimate goal 
of this research will be the higher performing inlets required for high speed flight. 
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National Aeronautics and 
Spoce Administration 

Lewis Research Center 


INTERNAL FLUID MECHANICS DIVISION NASA 


CFD FOR HIGH SPEED INLETS 




EXTERNAL COMPRESSION 



AXI-SYMMETRIC 


Notional Aeronoutics ond 
Spoce Administration 

Lewis Research Center 


INTERNAL FLUID MECHANICS DIVISION NASA 


THRUST DISTRIBUTION FOR SR-71 


MACH NO 

INLET 

ENGINE 

EJECTOR 

2.2 

13% 

73% 

14% 

3. + 

54% 

18% 

28% 
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Nationol Aeronautics and 
Space Administration 

Lewis Research Center 


INTERNAL FLUID MECHANICS DIVISION NASA 


CFD FOR HIGH SPEED INLETS 
1.3 < MACH <5.0 


PHENOMENA 

BUZZ 

UNSTART 

SHOCK/SHOCK INTERACTION 
OBLIQUE SHOCK/BOUNDARY LAYER 
NORMAL SHOCK/BOUNDARY LAYER 
GLANCING SHOCK/BOUNDARY LAYER 
CORNER FLOWS 


CODES 

3DPNS (PEPSIS) 

-BRILEY/MCDONALD LBI SCHEME 
3D-STEADY NS (PARC) 
-BEAMAVARMING ADI 
3D-UNSTEADY NS (PROTEUS) 
-BEAMAVARMING ADI 
-2nd ORDER TIME ACCURATE 


All of the above with 
BOUNDARY LAYER BLEED 
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CFD FOR HIGH SPEED INLETS 
OBLIQUE SHOCK BOUNDARY LAYER INTERACTION 



DENSITY GRADIENTS 
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Notionol Aeronautics and 
Space Administration 

Lewis Research Center 


INTERNAL FLUID MECHANICS DIVISION NASA 


SHOCK/BOUNDARY LAYER/BLEED INTERACTION 

EXPERIMENT/COMPUTATION COMPARISON 


o Experiment 


— Numerical results 


< 

> 

r°-- a ~° ° ° 

r- PPL/PT 

1 1 1 

/ 1 

1 1 I 


o 10 20 30 Ang=6.5 deg 

X [CM] 

WALL PRESSURE DISTRIBUTION 
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National Aeronautics and 

Space Administration INTERNAL FLUID MECHANICS DIVISION 

Lewis Research Center 


REFLECTED 



DISTANCE Y 


NASA 
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National Aeronautics and 
Spoce Administration 

Lewis Research Center 
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CFD FOR HIGH SPEED INLETS 

MACH 5.0 INLET 
STREAMLINE TRACES 
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Nationol Aeronautics and 
Space Administration 

Lewis Research Center 


INTERNAL FLUID MECHANICS DIVISION NASA 


CFD FOR HIGH SPEED INLETS 
SUMMARY 


Combined computational/experimental investigation has 
identified importance of 3D snock/boundary layer interactions 
for rectangular inlets at high speeds 


Possible control of glancing shock/boundary layer interaction 
through bleed now being explored. Little success to date. 


Understanding of three dimensional interactions may lead to 
better methods of flow control in inlets. 


ORIGINAL CONTAINS 
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THE USE OF A NAVIER- STOKES CODE IN THE WING DESIGN PROCESS 


S. Naomi McMillin 
NASA Langley Research Center 
Hampton, Virginia 23665-5225 

Abstract of Paper Presented at the 
NASA Computational Fluid Dynamics Conference 
NASA Ames Research Center 
March 6-9, 1989 


An ongoing investigation is being conducted in the Supersonic/Hypersonic 
Aerodynamics Branch at NASA Langley Research Center to determine the 
feasibility of incorporating the Navier- Stokes computational code, CFL3D, into 
the supersonic wing design process. The approach taken in this investigation 
is of two steps. 

The first step was to calibrate CFL3D against existing experimental data 
sets obtained on thin sharp-edged delta wings. The experimental data 
identified six flow types which are dependent on the similarity parameters of 
Mach number and angle of attack normal to the leading edge. The calibration 
showed CFL3D capable of simulating these various separated and attached- flow 
conditions . 

The second step was to use CFL3D to study the initial formation of 
leading-edge separation over delta wings at supersonic speeds. This study 
consisted of examining solutions obtained on a 65° delta wing at Mach number 
of 1.6 with varying cross-sectional shapes. Reynolds number was held constant 
at 1000000 and the Baldwin-Lomax turbulence model was used. The study showed 
that through the use of leading- edge radius and/or camber, the onset of 
leading- edge separation can be delayed to a higher angle of attack than 
observed on a flat sharp- edged wing. 

Based on the geometries studied, three wind-tunnel models are being 
designed to verify these results. These models are to be tested over a 
Reynolds number range of 2 x 10 b /foot to 8.5 x 10 b /foot. 
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The objective and approach for the present investigation. 

CFL3D code characteristics. 

Sketches of the wind-tunnel models tested by Miller and Wood in the 
Langley Unitary Plan Wind Tunnel (NASA TP -2430) . 

The types of flow classified from the Miller and Wood experimental 
test (NASA TP-2430). 

The computational test matrix at M * 2.8 superimposed on the chart 
which defines the flow types as functions of a and M normal to the 
leading edge. 

A comparison of computational results with experimental data for the 
75° delta wing at a = 8° and M - 2.80 (A1AA 87-2270). 

The major elements of the incipient separation study (second step of 
approach) . 

Computational results which quantify the effects of leading- edge 
radius at a = 4°. 

Computational results which quantify the effects of leading- edge 
radius at a = 8°. 

Computational results which quantify the effects of camber at a = 8°. 

A summary on the effects of leading-edge radius and camber on a 65° 
delta wing at M = 1.6. 

Major element of proposed wind-tunnel test. 

Concluding remarks. 
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OBJECTIVE 


To determine the feasibility of incorporating 
the Navier - Stokes computational code, 

CFL3D, into the supersonic wing design process. 

APPROACH 

Calibrate CFL3D with experimental data 
obtained on flat sharp-edged delta wings 
at supersonic speeds (AIAA 87-2270). 
Parameters: 

• leading edge sweep 

• angle of attack 

• Mach number 

Use CFL3D to study the initial formation of 
leading edge separation on delta wings. 
Geometric parameters: 

• leading edge radius 

• camber 
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CFL3D CODE CHARACTERISTICS 

• Developed at NASA-LaRC by Jim Thomas 

• Unsteady conservation-law form of compressible Navier-Stokes 
equations (thin layer laminar viscous model) 

• Laminar viscous modei or Baldwin-Lomax eddy viscosity 

turbulence model 

* 

• Upwind-biased spacial differencing implemented in finite 
volume form 

• Upwind (flux split) for convective terms and pressure 

• Central for shear stress/heat transfer 

• 2nd order accurate 

• Implicit time differencing 

• Streamwise relaxation/cross-flow ADI 

• 1st order accurate 

• Conical solutions obtained by 3-D solution 

• Single array of conically constructed volumes 

• Inflow equal outflow 

• Grids generated using an elliptic grid generation method, 

121 x 93 grid required to resolve details of flow patterns 


LEE-SIDE FLOW FLAT DELTA MODELS 
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FLOW CLASSIFICATIONS 



SHOCK-INDUCED 

SEPARATION 



NO SHOCK/NO SEPARATION 




SHOCK WITH NO 
SEPARATION 



SEPARATION BUBBLE 
WITH SHOCK 



SEPARATION BUBBLE 
WITH NO SHOCK 
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THEORETICAL TEST MATRIX 
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EFFECT OF VISCOUS MODEL ON NAVIER-STOKES SOLUTIONS 
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INCIPIENT SEPARATION STUDY 


OJ 



Examine effects of 

• Leading edge radius 

• Camber 

• Angle of attack 

• Reynolds number 

On the formation of 
leading edge separation 
on a conical 65° 
delta wing at M = 1.6 


Parametric cross-sections 


Leading edge radius: 

Sharp r/(b/2) = 0.0025 

(20:1 ellipse) 

e=— r — 

Camber (20:1 ellipse) 

4° camber 8° camber 

F~ = > E==— = 


r/(b/2) = 0.0100 

H 

10° camber 
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Crossflow Mach number: 



Surface pressure: 
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EFFECT OF LEADING EDGE RADIUS 

a = 8 r 


M = 1 .6 A= 65° Re = 1 000 000 
Turbulent boundary layer 
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I EFFECT OF CAMBER] 

EEH 

M = 1.6 A= 65° Re = 1 000 000 
Turbulent boundary layer 
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WIND TUNNEL MODELS 



Cross-sections: 

Sharp 

20:1 ellipse 

20:1 ellipse 
with 10° camber 





Wind Tunnel Test Description 


• M = 1 .6; a = 0° to 8° (initial conditions) 

• Surface pressure data (-100 ports) 

• Vapor screen photographs 

• Oil flows 


• Unitary Plan Wind Tunnel 

R e = 2 x 10 6 /foot, 5 x 10 6 /foot 

• 20" Supersonic Wind Tunnel 

R e = 2 x 10 6 /foot, 8.5 x 10 6 /foot 
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CLOSING REMARKS 


• The six distinct flow types observed over 
delta wings at supersonic speeds have 
been computationally reproduced using 
CFL3D (AIAA 87-2270) 


• CFL3D has been used to quantify the effects 
of leading edge radius and camber on leading 
edge separation. Wind tunnel test designed to 
confirm these results. 


N91 - 10 890 


APPLICATIONS OF A TRANSONIC WING DESIGN METHOD 


by 


Richard L. Campbell and Leigh A. Smith 
NASA Langley Research Center 


ABSTRACT 


The development of transonic computational fluid dynamic methods 
over the last fifteen years has had a tremendous impact on the design 
of airfoils and wings . A variety of approaches to automated design 
have been developed, such as the hodograph and fictitious gas methods 
for shock-free design, and inverse, numerical optimization, and 
predictor/corrector methods for the more general cases . Each of these 
techniques has advantages as well as limitations, and the designer 
must choose the one that best suits his specific task. 

A method for designing wings and airfoils at transonic speeds 
using a predictor/corrector approach has been developed. The pro- 
cedure iterates between an aerodynamic code, which predicts the flow 
about a given geometry, and the design module, which compares the cal- 
culated and target pressure distributions and modifies the geometry 
using an algorithm that relates differences in pressure to a change in 
surface curvature . The modular nature of the design method makes it 
relatively simple to couple it to any analysis method. The iterative 
approach allows the design process and aerodynamic analysis to con- 
verge in parallel, significantly reducing the time required to reach a 
final design. Viscous and static aeroelastic effects can also be 
accounted for during the design or as a post-design correction. 

Results from several pilot design codes indicated that the method 
accurately reproduced pressure distributions as well as the coordi- 
nates of a given airfoil or wing by modifying an initial contour. The 
codes were applied to supercritical as well as conventional airfoils, 
forward- and aft -swept transport wings, and moderate-to-highly swept 
fighter wings . The design method was found to be robust and effi- 
cient, even for cases having fairly strong shocks. Comments from a 
user in industry indicated that for a specific design problem, this 
design method was about 25 times faster than a numerical optimization 
approach that utilized the same aerodynamic analysis code. 
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Flow conditions 
Initial geometry 


Aerodynamic 
analysis module 


Aeroelastic 
analysis module 


Pressure distribution 
current geometry 
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HYBRID DESIGN ALGORITHM 


• Subsonic regions 

A C =A C p A (1 +C 2 ) B 


• Supersonic regions (M>1.15) 


AC = °' 05A i (ACp) 

( 1 +( y ') 2 ) 1 ' 5 dx 
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SUPERCRITICAL AIRFOIL DESIGN 


M =0. 77 


a= 1 .0 




initial 

FINAL 

♦♦♦♦♦♦ TARGET 
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SUPERCRITICAL AIRFOIL DESIGN 



— INITIAL 

— FINAL 
—TARGET 


x/c 
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NYU CODE- PARTIAL AIRFOIL MODIFICATION 



Inital 

Final 

+++++++ Target 


x/c 
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NYU CODE- PARTIAL AIRFOIL MODIFICATION 



0 .2 .4 .6 .8 1.0 


x/c 
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EXECUTIVE JET DESIGN 
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EXECUTIVE JET DESIGN 


M = 0.78 
T| = .12 



— Analysis 

— Final 

++++++ j ar g e t 
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EXECUTIVE JET DESIGN 

M = 0.78 
T| = .12 


. 08 1 - 



0 .2 .4 .6 .8 

x/c 


— Final 


.0 




HIGHLY SWEPT FIGHTER 


WING GEOMETRY 
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HIGHLY-SWEPT FIGHTER WING CASE 


c 


p 


M = 0.85 




— Initial 

— Final 
*♦♦♦♦♦ Target 


0 .2 .4 .6 .8 1.0 

x/c 


HIGHLY-SWEPT FIGHTER WING CASE 

M = 0.85 



— Initial 

— Final 
♦♦♦♦♦♦ Target 
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HIGHLY-SWEPT FIGHTER WING CASE 

r| = 0.77 



Initial 

Final 

Target 
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HIGHLY-SWEPT FIGHTER WING CASE 

T1 = 0.52 



Initial 

Final 

Target 


x/c 
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CONCLUDING REMARKS 


•A predictor/corrector method for designing airfoils and 
wings at transonic speeds has been developed 

•Hybrid algorithm used in strongly supercritical regions 
to improve convergence of design process 

•Predictor/corrector approach allows geometry 
constraints to be easily included 

•Results from airfoil and wing design pilot codes 
indicate that the method is robust, efficient, and 
accurate for a wide variety of configurations 

•Techniques for accounting for viscous and aeroelastic 
effects have been included 
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ALGORITHMS AND TOOLS 
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Joseph L. Steger 
Senior Staff Scientist 
Fluid Dynamics Division 
NASA Ames Research Center 
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An Embedded Grid Formulation 
Applied to Delta Wings 


James L. Thomas 

NASA Langley Research Center, Hampton, Virginia 
Sherrie L. Krist 

Vigyan Research Associates, Inc., Hampton, Virginia 


An embedded-grid algorithm tor the Euler/Navier-Stokes equations is developed and applied to delta wings at high angles of attack in low- 
speed flow. The Navier-Stokes code is an implicit, finite-volume algorithm, using flux-difference-splitting for the convective and pressure 
terms and central differencing for the viscous and heat transfer terms. Calculations are compared with detailed experimental results over an 
angle of attack range up to and beyond the maximum lift coefficient, corresponding to vortex breakdown at the trailing edge, for a delta 
wing nominally of unit aspect ratio. The results indicate that the overall flowfield, including surface pressures, surface streamlines, and 
vortex trajectories, can be simulated accurately with the global grid version of the present algorithm. However, comparison of computed 
velocities and vorticity with experimentally measured off-body values at an angle of attack of 20.5 degree indicates the core region is 
substantially more diffuse in the computations than that measured with either a five-hole probe or a laser velocimeter. 

Embedded grids, used in the present work to improve the numerical discretization in the core region, are formulated within the framework 
of the implicit, upwind-biased multi-grid algorithm. Structured levels of local nested refinements are made. The refinements, which should 
be based on the local truncation error of the solution, are currently based on the location of the vortex core trajectory. The embedded grid 
boundary data are determined through linear interpolation of the coarser underlying grid. The values of flow variables in the coarser grid 
cells underlying finer grids are determined by a volume-weighted restriction of the finer grid values. To compensate for the limited extent 
of the implicit operators in the embedded grids, each level of grid refinement uses corrections from coarser grids, which can be either 
embedded or global grids, to accelerate the convergence to the steady state of the embedded blocks. 

Three-dimensional results for both Euler and Navier-Stokes calculations are shown, with up to 3 levels of embedded refinement. The 
embedding procedure was effective in eliminating a crossflow secondary separation produced in the Euler solutions on coarse grids. The 
core velocities in the Euler solutions are substantially higher than in the viscous solutions and in both, the embedded grids can efficiently 
reduce the numerical diffusion in the flowfield. The computed core velocities in the viscous calculations are increased on the embedded 
grids but are lower than the measured values with the grids used in the present study. 
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INTRODUCTION 


Motivation : 

Performance and stability of aircraft at high angles of attack 
are dominated by the nonlinear interaction of shed vortices 


Present objective : 

Compare predictions of Navier- Stokes code with detailed 
experimental data for delta wings with both global and 
embedded grid schemes 


Primary vortices 


Trailing edge 


Secondary vortex 



Apex 
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OUTLINE OF PRESENTATION 


• Introduction 

• Global grid scheme 

- Computational algorithm 

- Comparisons with experiment 

• Embedded grid scheme 

- Spatial differencing 

- Time advancement 

• Computational results 

- Euler 

- Navier- Stokes 

• Concluding remarks 
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GLOBAL GRID ALGORITHM 

Upwind-biased Navier-Stokes code CFL3D 


dQ | dF | dG | d{H^-Hy) 


dt d£ dr] d( 


• Time-dependent conservation law form of compressible Navier- 
Stokes equations 

• Implicit (diagonalized) time advancement algorithm 

• Upwind-biased spatial differencing 

o Flux difference splitting for convective/pressure terms 
o Central differencing for shear stress/heat transfer terms 

• Thin-layer with algebraic turbulence model 

• Semi-discrete finite- volume implementation 

• Zonal (patched) grids 
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COMPARISONS WITH EXPERIMENTS 

Delta wings at high alpha in low speed flow 


• Hummel (DFVLR) AR=1 wing 

- R l = 0.95 x 10 6 

- Force and moment 

- Surface oil flows 

• Kjelgaard and Sellers (BART) A = 75° wing 

- R l = 0.5 - 1.5 x 10 6 

- Surface oil flows 

- Flowfield measurements ( 5-hole probe / LDV ) 
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SECONDARY AND TERTIARY SEPARATION LINES 
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AR = 1 DELTA WING 
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Embedded Grid Scheme 



Coarse grid 

Fine grid 

• Volume-weighted restriction of 
finer grid to coarser grid cells 

• Linear interpolation for projections 
of finer grid cells into coarse grid 
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Embedded Grid Scheme 


Level 


Embedded Grid 


Global Grid 


Global Grid 



E: Navier-Stokes Calculation 
R: Restriction of Residual and Q 
P: Prolongation of Correction 
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PARTICLE TRACES 


33 x 33 x 19 Mesh 



w/Grid Embedding 
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EULER SOLUTION FOR 76 DEGREE DELTA WING 

u velocity component 



/ 
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GRID REFINEMENT EFFECT 
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3-D NAVIER STOKES RESULTS 


a = 20.5 °; 

A = 75 °; R l = 

0.5 x 10 6 

Mesh 

Grid points 

u max / u oo 



@ x/L = 0.7 

33x33x19 global 

20,691 

1.20 

33x33x19 global 

233,160 

1.80 

w/ 3 embedded levels 



65x65x37 global 

156,325 

1.55 

65x65x37 global 

512,770 

1.92 


w/ 1 embedded level 
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STREAMWISE VELOCITY CONTOURS 

x/L = 0.7 ; a = 20.5 


Experiment 
Umax/voo ~ 2.25 


Navier- Stokes 




377 


CONCLUDING REMARKS 


• Global grid algorithm predicts overall flow topology 
and surface characteristics of delta wings in good 
agreement with experiment 

- Force and moment 

- Surface oil flows 

- Vortex trajectories 

• Embedded grid algorithm provides efficient way to predict 
local off-surface flows 


- Vortex core velocities 

- Downstream interactions 
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Unstructured Mesh Solution 
of the Euler and Navier-Stokes Equations 


Timothy J. Barth 
CFD Branch 

NASA Ames Research Center 
Moffett Field, CA 94035 


Mesh generation procedures as well as solution algorithms for solving 
the Euler and Navier-Stokes equations on unstructured meshes are pre- 
sented. The solution algorithms discussed utilize approximate Riemann 
solver, upwind differencing to achieve high spatial accuracy. Numerical 
results for Eider flow over single and multi-element airfoils are presented. 
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Motivation/Requirements 


• Generality 

• Mesh Generation of Complex Geometries 

• Inviscid/ Viscous, Subsonic/Transonic Flows 

• Accuracy 

• High Spatial Accuracy for Smooth Flows 

• Mesh Adaptation 

• Shock Capturing for Discontinuous Flows 

• Generalized Upwind Scheme 

• Efficiency 

• Edge Data Structure 


• Edge and Face “Coloring” 


Mesh Generation 


Structured 

• Algebraic, Elliptic, Hyperbolic, Parabolic, 
Delaunay Triangulation 

• 0[n log n) 

• Maximize Minimum Angle 
Advancing Front 

• Probably 0(n logn) 
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Fig. 2.1a) Three-component configuration. 


Fig. 2.1b) Grid about main element. 



Fig. 2.1c) Grids plotted atop one another. 


Fig. 2. Id) Grid after elimination of unwanted pts. 



Fig. 2.1e) Grid after reconnection of points. 


Fig. 2. If) Grid after reconnection, detail. 


Figure (2.1) Mesh generation synthesis of 3 element airfoil. 
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Finite- Volume Spatial Scheme 


Euler equations in integral form: 

d_ 
dt 


I u 


da + (p f(n) dl 
n J dci 


= 0 


where f(n) = n x f + n y g 


Numerical approximation about polygonal face: 

e r d{fj) 

— u da=~Y] Al t h(u+,u“ ; n,), j 

dt Ja fj i=1 


Roe Flux Function: 


h( u+ , u ; n) =-(f(u + ; n) + f(u ;n)) 

- ^|A( u+ , u “;n)|(u + - i 


1 ,n{f) 


~) 
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Reconstruction - Evolution 


• Solution Reconstructions 

• Piecewise Constant 

• Piecewise Linear u(x,y) = u(xo,yo) + 3>Vu • A r 

• Gradient Estimation in each face, Vu = — f ao un dl 

’ aft Jail 

• Monotonicity Enforcement 

• Limiting (Multidimensional), G [0, 1] 

\ 

• Solution Evolution 

• Three Stage Runge-Kutta Time-Stepping 
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Numerical Examples 


• NACA 0012 airfoil = .8, a = 1.25°) 

• Subdivided Structured Mesh (193x33) 

• Irregular Mesh 

• 3 Element Airfoil = .2, a = 0.0°) 
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3 -Element Airfoil Cp Distribution 



2 - 
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.2 


x location 
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Current and Future Directions 


• Quadratic Reconstructions 

• Viscous Calculations with Turbulence Modelling 

• Adaptive Mesh Refinement 
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3-D UNSTRUCTURED GRIDS 
FOR THE SOLUTION OF 
THE EULER EQUATIONS 


C. Gumbert 

NASA Langley Research Center 

P. Parikh, S. Pirzadeh 

Vigyan Research Associates 

R. Lohner 

George Washington University 


NASA CFD Conference 
NASA Ames Research Center 
March 7-9, 1989 
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Advancing Front Grid Generator 

The advancing front technique is being used to develop a code 
to generate grids around complex three dimensional configurations for 
use in computing the invisid flow solutions by the Euler equations. By 
the advancing front technique points are introduced concurrently with 
the connectivity information so that a separate library is not required. 
The generation of a 3-D grid is accomplished in several steps. First 
the boundaries of the domain to be gridded must be described by 
two-, three— or four-sided surface patches. Next, a background mesh 
is required to control the grid spacing and stretching throughout the 
domain. This coarse tetrahedral grid is not required to conform to 
any of the boundaries. Next, each of the patches is mapped to 2-D, 
triangulated by the advancing front technique and mapped back to 
3-D. These triangles form the initial front for the generation of the 
final tetrahedral mesh. 
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ADVANCING FRONT GRID GENERATOR 


• The field points are introduced while grid is generated. 

- Does not require a separate library to introduce points before 

grid generation. 

- Allows adaptive regeneration of grids easily. 

• Generation Steps: 

- Define boundaries of the domain using surface patches. 

- Set up background grid to control local grid characteristics. 

- 2-D tri angulation of patches 
=$► Initial Front 

- Advancing front in 3-D to fill region with tetrahedra. 


Advancing Front in 2-D 


The figure is a schematic showing the advancing front technique 
in 2-D. The initial front is a set of line segments called faces. In 3-D 
a face is a triangle. In succession, each face of the front is deleted 
and a new point in the field is introduced. As in the center figure, if 
the new point is close to an existing point, then the existing one is 
used instead. The front advances until it closes in on itself at which 
point the region is fully gridded. 
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ADVANCING FRONT IN 2-D 



initial front 



Surface Grid for F-18 Configuration 


The figure shows the surface grid for an F-18 fighter configuration. 
The grid was generated for only one half of the configuration but 
it has been mirrored in the picture. The grid consists of 367,000 
tetrahedrons using nearly 66,000 points of which 10,000 lie on the 75 
surface patches defining the airplane and the computational box. It 
required about 400 seconds on the NAS CRAY-2 computer. 
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Finite Element Flow Solver 

The flow solver developed in conjunction with the grid generator 
uses a two-step Taylor Galerkin finite element method for the Eu- 
ler equations cast in Arbitrary Lagrangian-Eulerian (ALE) form. The 
Galerkin weighted residual method is used to perform the spatial dis- 
cretization. Timestepping options are available for steady-state or tran- 
sient problems. Accurate solutions without spurious over/undershoots 
can be obtained using Flux Corrected Transport (FCT) techniques. 
Second order pressure or Lapidus damping is used near shocks. Adap- 
tive mesh refinement is used to better capture sharp gradients. This 
technique has not been fully implemented in 3-D. 
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FINITE ELEMENT FLOW SOLVER 


• Two-step Taylor-Galerkin scheme 

• Spatial discretization performed 

via Galerkin weighted residual 
method using linear elements 

• Options for: 

Global/local timestepping 

Second order pressure or Lapidus damping 

Flux Corrected Transport (FCT) 

Adaptive H-reflnement 
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Results In 2-D 

As a means of determining the accuracy of this code, results 
are presented for two 2-D flow calculations. The first case is the 
transonic flow calculation for an NACA 0012 airfoil at = .80 
and a = 1.25°. Comparison is made with results from FL052, a 
finite volume method calculated on an O-type grid. 

The second case is the flow field around a 20° ramp at = 3.0. 
The exact solution is known and can be compared. This case also 
illustrates the use of adaptive mesh refinement. 
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RESULTS IN 2-D 


• Comparison of NACA 0012 airfoil 

at Moo = -80, a = 1.25° 

FEFL027 ( unstructured mesh, Finite 
Element solver) 

FL052 (O-grid, finite volume method) 

• 20° ramp flow 

at Moo = 3.0 

adaptive mesh refinement 
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Grid for NACA 0012 Airfoil 

The figure shows the unstructured grid used for the calculation. 
There are nearly 5000 points in the grid of which 266 points lie on 
the airfoil surface. This is similar to the 193 points in the structured 
grid of FL052. 
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Pressure Coefficient Comparison 
for NACA 0012 Airfoil 

The pressure coefficient distributions are compared in the figure. 
The results of the calculation made on the structured grid are des- 
ignated by diamonds connected by straight lines. The results of the 
calculation on the unstructured grid are designated by filled circles. 
The solutions are very similar. The location and strength of the nor- 
mal shock on the upper surface are nearly identical. The finite ele- 
ment solver captured the shock more sharply than did FL052. The 
finite element solution shown was obtained in 4000 steps which used 
1.9 x 10 -3 seconds/point/timestep on a Convex C2 computer. 
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Exact Solution For Supersonic 
Flow Over A Ramp 

This figure shows a schematic of the second case, supersonic flow 
over a 20° wedge. The exact solution is shown for the regions 
designated 1, 2, and 3. 
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EXACT SOLUTION FOR 
SUPERSONIC FLOW OVER RAMP 



e =20° 


Grid Adaptation for 20° Ramp 

This figure shows the final adapted grid for the solution. The levels 
of adaptation are shown by different colors. 
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Mach Number and Pressure 
Contours for 20° Ramp 

The calculated Mach number and pressure are shown in the fol- 
lowing figures. The shock angle and the conditions in the regions 
designated 1, 2, and 3 are predicted accurately. The figure also shows 
the interaction of the oblique shock with the expansion fan and the 
reflection of the shock off the straight wall. 
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Surface Pressure 
Contours on F-18 Aircraft 

The figure shows the pressure contours on the surface of an F- 
18 aircraft configuration at Moo = 0.3 and a = 15°. The solution 
was calculated on a grid for which the engine inlet and nozzle were 
blocked. This solution was run 600 iterations in 159 minutes on the 
NAS CRAY-2 computer to decrease the maximum residual by two 
orders of magnitude. 
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Pressure Contours For Shock 
Emanating From A Corner 

Sharp gradients in density and pressure were imposed at the in- 
tersection of three perpendicular walls as the initial conditions for a 
transient flow calculation. The figures show the pressure contours on 
the walls at two different time steps. The grid on the three walls is 
also shown. Additional grid points were added during the calculation 
to better capture the sharp gradients. No derefinement was done after 
the large gradients had passed. 
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Future Plans 

Further development of the grid generator will tend along three 
paths: (1) implementation of existing CAD sytem for faster and easier 
surface definition, (2) enhancement of the surface element library to 
allow slope continuity across surface patch boundaries and (3) faster 
and better ways of defining the background grid. 

Further development of the flow solver is expected to involve 
incorporation of the viscous terms into the flow equations, use of 
implicit or semi-implicit algoriths and implementation on parallel 
computers. 
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FUTURE PLANS 

• For the grid generator: 

- Link to commercial CAD/CAM systems for faster and 
more accurate surface definition. 

- Faster and better ways to input background grid information. 

- Enhancement of surface element library. 

• For the flow solver: 

- Incorporation of viscous terms. 

- Implicit or semi-implicit flow solver for transient problems. 

- Implementation on parallel computers. 
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FLUX SPLITTING ALGORITHMS FOR TWO-DIMENSIONAL 
VISCOUS FLOWS WITH FINITE-RATE CHEMISTRY 


Jian-Shun Shuen 

Sverdrup Technology , Inc. 

NASA Lewis Research Center, Cleveland, Ohio 4J185 


Meng-Sing Liou 

NASA Lewis Research Center, Cleveland, Ohio 44 185 


ABSTRACT 

The Roe flux-difference splitting method has been extended to treat two-dimensional viscous 
flows with nonequilibrium chemistry. The derivations have avoided unnecessary assumptions or 
approximations. For spatial discretization, the second-order Roe upwind differencing is used for 
the convective terms and central differencing for the viscous terms. An upwind-based TVD scheme 
is applied to eliminate oscillations and obtain a sharp representation of discontinuities. A two-stage 
Runge-Kutta method is used to time integrate the discretized Navier-Stokes and species transport 
equations for the asymptotic steady solutions. The present method is then applied to two types of 
flows: the shock wave/boundary layer interaction problems and the jet in cross flows. 


N91-10894 


438 


OBJECTIVES 


1. To develop a numerical algorithm for real gases which is robust for high 
Mach number flows and has the ability to accurately capture strong shocks 
and other types of discontinuities. 

2. To include nonequilibrium chemistry in the formulation. 

3. To make the formulation consistent and simple. 
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GOVERNING EQUATIONS 


The 2-D, thin layer Navier-Stokes and species equations for a chemically react- 
ing gas of N species: 


au 3F(u) aG(u) aw(u,Uy) 

dt + dx + dy dy -S 

where 
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THERMODYNAMIC AND TRANSPORT MODELS 


Equation of State 


P 


P[.Pi ^1 > ^2> ••••(-' ff_ x ) 



Physical Properties 

• Mass diffusion is approximated by Fick’s law. Binary mass diffusivity ob- 
tained from Chapman-Enskog theory. 

• For the individual species: C p , k, and /x are determined by fourth-order 
polynomials of temperature. 

• For the gas mixture: C p obtained by mass concentration weighting of each 
species, k and /x obtained using Wilke’s mixing rule. 

Chemistry Model of Air 

• Five (5) species (O 2 , N 2 , O, N, NO) and eleven (11) elementary reaction 
steps. 

• Neglecting ionization as well as thermal nonequilibrium. 
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SOLUTION METHOD 


1. Two-stage Runge-Kutta method for time integration of the governing equa- 
tions. 

2. Chemical source terms treated implicitly so that the stiffness associated 
with the source terms will not upset numerical stability. 

3. Second-order upwind differencing for convective terms, central differencing 
for viscous terms. 

4. Upwind-based TVD scheme to eliminate oscillations and to obtain a sharp 
representation of discontinuities. 
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UPWIND DIFFERENCING 


• Good Property for Capturing Discontinuities. 

• Flux- Vector Splittings 

1/ Steger- Warming Splitting 
1/ Van Leer Splitting 

• Flux-Difference Splittings 

/ Roe Splitting 
Osher Splitting 
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ROE FLUX-DIFFERENCE SPLITTING 


To construct Roe flux-difference splitting, one usually defines an average state 

A 

U such that 

AF = AAU, AG = BAU, A(-) = (•)„- ( ) t , 

A = A(U), B = B(U), U = U(U L ,U B ) 

The eigenvalues of the Jacobian A are (for N=5): 

A(A) = u, u, u — a, u + a, u, u, u, u 

A is diagonizable 

A = S a A a S“ 1 , diag A a = A(A). 

Splitting of the diagonal matrix A a results in the splitting of A^: 

A = A++A-, A ± = s x Ajs ; 1 

Thus, the splitting of the flux difference AF: 

AF = AF + + AF“, AF ± = A ± (U)AU. 

A A / 

Goal: Find the average state U = U(Ul,Ur). 
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ROE’S AVERAGE STATE 


• AF = A(U)AU 

• For a calorically perfect gas: the Roe average state is easily obtained. 

• For a real gas: the simplicity is lost. 

Choose the set of average quantities 

(p, u, u, e, , Y{ ,t=l,4> PpiPei Pc- • . . ) 

1,1 — 1,4 

Define the Roe-average operator p, 

/ = M(/) = ^ 7 7^- . d=(p L /p R )*, f = u,v,e,H,Yi, P = P R d. 

1 + a 

The condition left to be satisfied is: 

N - 1 

A p = ppAp + p e Ae + Y2 Pc ( A C{. (A) 

t=i 

Condition (A) contains N+l variables for only one equation. 


445 


ROE’S AVERAGE STATE (cont’d) 


Q — (Pe, Pp, p Gi ), Pp — Pp(p, £, C»), etc. 
Q = {Pe,Pp,P c .) 





(Pe, Pp> P c .) = (p e , Pp, P c .) for a calorically perfect gas. 
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NUMERICAL TEST 


• Constituent elements of the present flux-splitting method validated by com- 
parison with exact solutions or experimental data in our previous work on 
1-D Euler equations for reacting gases, and 2-D Navier-Stokes equations 
for ideal gases. 

• Direct comparison of the present 2-D reacting flow N-S calculations with 
experimental data is hindered by the lack of suitable data for high Mach 
number reacting flows. 


Two classes of problems are considered in the following: 

(1) shock wave/laminar boundary layer interaction problems. 

(2) jet in cross flow problems. 
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SHOCK WAVE/BOUNDARY LAYER INTERACTION 
WALL PRESSURE AND TEMPERATURE 

Moo = 8.0, 0 = 22°, R e = 1.57 x 10 6 



o 



JET-IN-CROSS FLOW 


VELOCITY VECTORS AND WALL PRESSURE 

Moo = 6.0, R e = 3.21 x 10 6 
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DUAL JETS-IN-CROSS FLOW 


VELOCITY VECTORS AND WALL PRESSURE 
Moo = 6.0, R t = 3.21 x 10 6 
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SUMMARY 


1. New average state of Roe’s splitting for chemical reacting gases is proposed. 

2. A numerical method based on Roe flux-splitting for solving the Navier- 
Stokes and species equations for chemical nonequilibrium gases has been 
developed. 

3. Numerical results show the present method has the potential for calculating 
high Mach number flows with nonequilibrium chemistry and complex wave 
structure. 

4. Implementation of an implicit time solver is in progress. 

5. Thermal nonequilibrium is being considered in the next stage of code de- 
velopment. 

6. Work for further code validation is needed. 
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Visualization of Fluid Dynamics at NASA Ames 

Val Watson 

NASA Ames Research Center 


The hardware and software currently used for visualization of fluid dynamics at 
NASA Ames is described. The software includes programs to create scenes (for 
example particle traces representing the flow over an aircraft), programs to 
interactively view the scenes, and programs to control the creation of video tapes 
and 16mm movies. The hardware includes high performance graphics 
workstations, a high speed network, digital video equipment, and film recorders. 

With the current workstations, a scientist can interactively view flow over 
simplified objects, such as the flow over a circular cylinder. For complex objects, 
such as an aircraft, the workstation creates each picture too slowly to gain a sense of 
the dynamics of the flow. Therefore, each picture is stored frame by frame on a 
video tape or 16mm film and then the video or movie is played back at normal 
speed to illustrate the flow dynamics. 

The upgrade in workstations planned for this year is expected to permit moderately 
complex pictures (pictures that can be represented by 10,000 polygons or less) to be 
created at a rate of 10 frames per second — fast enough to gain a sense of the flow 
dynamics. Therefore, these workstations should permit interactive viewing of the 
flow over complete aircraft rather than just simple objects. Upgrades planned this 
year for software should provide a more effective interface for controlling the 
interactive viewing. 

A comparison of the upgrades planned this year with an ideal simulation and 
visualization environment shows that there is still potential for major 
improvements in both software and hardware. The greatest potential for 
improving the environment is the development of software to extract and illustrate 
the essence of very complex phenomena. 

Results presented by other scientists during this conference clearly demonstrate the 
effectiveness of the current visualization tools for assisting in the understanding of 
complex simulations, but it is also clear that we are a long way from utilizing 
visualization tools to their full extent. 
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Outline 


NASA’s current visualization tools 
Capabilities with the current tools 
Upgrades planned for this year 
Potential for further improvements 
Conclusions 
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Visualization Hardware 
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Video Recording System 
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Software 


Simulation 

Visualization and Recording 


Scene 

Creation 

Scene 
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Animation 
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Current Capabilites in Visualization 


Interactive viewing with workstations 

Dynamic illustration of wire frame objects 

Dynamic illustration of simple solid objects 

Static illustration of complex solid objects 

Playback viewing with video or film 

Dynamic illustration of complex solid objects 
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Upgrades in Workstations Planned this Year 


Basic Workstation features 
Central processor from 2600 to 2000 Dhrystones 
Arith processor from 0.1 to 6.0 MFLOPS 
Primary memory from 4 to 32 MBytes 
Secondary mem. from 474 to 1000 MBytes 
Graphics Features 

3D coord, transf. from 80K to 800K coord./sec 
Solids rendering from 0.5K to 100K polygons/sec 

( = 10,000 polygons at 10 frames/sec) 
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U pgrades in Software Planned this Year 

Combination of old modules into a single program 
Changes to take advantage of new workstations 

Changes to make visualization more interactive 
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Comparison of Current Capabilities with "Ideal" 

for Vision 


Feature 

Current 

"Ideal" 

spatial res. 

1 minute of arc 

of the same order 

§ color res. 

1 6 million colors 

of the same order 

freq. resp. 

(scene dependent) 

15 frames/sec 

field of view 

1/5 steradian 

5 steradians 
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Comparison of the Ideal vs Current Field of View 


Field of view of the human eyes 

Field of view of 19" display 





Comparison of Current Capabilities with "Ideal" 

for Interactive Control 


Current 

mouse and keyboard 
"Ideal" 

6-degree of freedom control 
voice recognition 
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Potential for Improvements in Software 

Tools for extraction of critical features 

Tools for highlighting critical features and 
suppressing less important features 
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CONCLUSIONS 


The NAS decision to invest in visualization tools 
has been justified by improved analysis capabilities 

The most effective use of vis. tools is the routine 
interactive viewing of solutions and techniques 

Current workstations permit interactive viewing of 
the flows over simple shapes 

Workstations to be procurred should permit 
interactive viewing of flows over complex shapes 
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CONCLUSIONS (continued) 


A new generation of software is being developed 
to take advantage of new workstation capabilties 
to make the visual analysis more interactive 

The most critical task is overcoming visual clutter 
developing techniques to extract the "essence" 
developing techniques to illustrate the "essence" 


We are a long, wav from utilizins our present 
visualization tools completely 


465 


For Further Information 


NAS procurement for the new workstations 

Rosemary Buchanan 
Mail Stop 258-6 
NASA Ames Research Center 
Moffett field, CA 94035 
(415) 694-4610 

Scientific visualization materials 

Stanford report on automated feature extraction 
1988 CFD Highlights video 
Stereo slides 

Pat Mumford Elson 
Mail Stop 258-2 
NASA Ames Reseach Center 
Moffett Field, CA 94035 
(415) 694-4463 
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COMPUTATIONAL FLUID DYNAMICS ON A 
MASSIVELY PARALLEL COMPUTER 

Dennis Jespersen, CFD Branch, NASA/Ames 
and 

Creon Levit, NAS Applied Research Branch, NASA/Ames 


Abstract. We have implemented a finite-difference code for the compressible 
Navier- Stokes equations on the Connection Machine, a massively parallel computer. 
The code is based on the ARC2D/ARC3D program and uses the implicit factored 
algorithm of Beam and Warming. The code uses odd-even elimination to solve 
linear systems. We give timings and computation rates for the code, and compare 
with a Cray XMP. 
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COMPUTATIONAL FLUID DYNAMICS ON A 
MASSIVELY PARALLEL COMPUTER 

Dennis Jespersen, CFD Branch NAS A/ Ames 

and 

Creon Levit, NAS Applied Research Branch, NASA/Ames 


• Goal: implement a real CFD code (ARC2D/ARC3D) on 
a massively parallel computer, the Connection Machine. 
> Implicit 


Is Connection Machine useful for CFD? 
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OUTLINE 


• CM architecture 

• ARC2D 

• *Lisp 

• Odd-even elimination 

• Current status and results 

• Problems and prospects 
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CM ARCHITECTURE 


I 


• Very large number of physical processors (up to 64K) 
« SIMD 

• Single-bit processors 

• 2K words (32 bits each) memory per processor 

• Floating point accelerator available 

• Context bit (each processor on or off) 

• Peak performance (computational kernels) 2-4 Gflops 
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CONNECTIONS 


® 16 processors per chip 

• Chips wired in 12-dimensional hypercube 

> If distance=power of 2, communication in 2 hops 
@ General router 

® Nearest-neighbor (NEWS) available and fast 

> (Reminiscent of ILLIAC) 

• Wrap-around at edges 
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VIRTUAL PROCESSORS 


® Each physical processor can be configured as 1, 2, 4, 8, . . . 
virtual processors. 


VP ratio 1 2 4 8 16 32 64 

memory /processor (words) 2K IK 512 256 128 64 32 


• Transparent to user. 

® Arithmetic slows down (slightly slower than linear) 

• Communication relatively faster 

> Memory references that were off-chip may become 
on-chip. 


ARC2D 


2D Compressible Navier- Stokes 

d±Q + d^E + d^F — Re + d^Fy) 

Central space differencing, implicit Euler in time, fac- 
tored, delta form: 

(I + A t8^A n - A t~8^R n )(I + A t8 v B n - At^-8 v S n )AQ 

Re Re 

= - A t(8^E n + 8 v F n - + 6 V F?)) 

Pulliam-Chaussee diagonal form: 

T^(I + AtS^A^ - At—8^A^ V ) 

N(I + AtSrjArj - At^—6 ri A ri<v )Tr 1 AQ n = rhs 

Re 5 ^ 
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ARC2D on CM 


• One grid point per (virtual) processor. 

• Write in *Lisp 

> No line-by-line translation 

• Shear-layer 

• Cylinder 

• Typical size 128 x 64 (8K processors, VP ratio = 1) 

• CM well-suited for explicit time-step methods (e.g., Runge- 
Kutta) 


IMPLICIT SOLVER: ODD-EVEN ELIMINATION 


Closely related to cyclic reduction 
One row corresponds to one processor 
Communication is always to processors 2 k distant 
2 n x 2 n system solvable in n stages 
Algorithm makes no sense on uniprocessor machines 
Stability not a problem in practice 
> Some stability results provable, e.g. B(— At, 1, At) 


476 


ODD-EVEN ELIMINATION: ALGORITHM 
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ODD-EVEN ELIMINATION:PERIODIC 


• Same code 

• Same cost 

• Periodic solves no more costly than nonperiodic solves 
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PROGRESS 


• 2D Explicit (Runge-Kutta) 

• 2D Implicit (scalar and block tridiagonal) 

> Shear layer 

> Cylinder (low Re) 

• 3D Explicit (Runge-Kutta) 
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TIMINGS 

• Timings, Mflop rates (16k processors, implicit scalar tridi- 
agonal algorithm) 


Grid size 

VP ratio 

Mflops 

sec/step XMP time/CM time 

128 x 128 

1 

58 

0.43 

0.66 

256 x 128 

2 

86 

0.60 

0.75 

512 x 512 

16 

106 

4.2 

0.86 

Timings, Mflop rates 

(16k processors, explicit algorithm) 

Grid size 

VP ratio 

Mflops 

sec/step XMP time/CM time 

128 x 128 

1 

63 

0.43 

0.69 

256 x 128 

2 

123 

0.44 

1.34 

512 x 512 

16 

241 

1.8 

2.62 
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TIMINGS: 2D BLOCK IMPLICIT AND 3D 

2D Timings, 16k processors, block tridiagonal algorithm 


Grid size VP ratio sec/step XMP time/CM time 
128 x 128 1 2.61 0.23 

256 x 128 2 3.21 0.37 

256 x 256 4 4.84 0.49 

512 x 256 8 9.66 0.49 

• 3D Timings, 8k processors, explicit algorithm 

Grid size VP ratio sec/step 
32 x 32 x 16 2 1.94 

64 x 32 x 32 8 3.80 

64 x 64 x 32 16 6.53 


PROBLEMS AND PROSPECTS 


Memory limitation 
Communication bottleneck? 

Boundary conditions 

> Complicated boundary conditions - inefficient? 
Future plans: 

> 3D Implicit 


SESSION XII 


HYPERSONICS / AFE 


Chairman: 

Chien P. Li 

Advanced Programs Office 
NASA Johnson Space Center 
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Conservation Equations and Physical Models for Hypersonic Air 
Flows Over the Aeroassist Flight Experiment Vehicle 

Peter A. Gnoffo 

NASA Langley Research Center 
Hampton, VA 23665 
Mail Stop 366 


Abstract 

The code development and application program for the Langley Aerothermodynamic 
Upwind Relaxation Algorithm (LAURA), with emphasis directed toward support of the 
Aeroassist Flight Experiment (AFE) in the near term and Aeroassisted Space Transfer Vehicle 
(ASTV) design in the long term is reviewed. LAURA is an upwind-biased, point- implicit 
relaxation algorithm for obtaining the numerical solution to the governing equations for three- 
dimensional, viscous, hypersonic flows in chemical and thermal nonequilibrium. The algorithm 
is derived using a finite-volume formulation in which the in viscid components of flux across 
cell walls are described with Roe’s averaging and Harten’s entropy fix with second-order correc- 
tions based on Yee’s Symmetric Total Variation Diminishing scheme. Because of the point- 
implicit relaxation strategy, the algorithm remains stable at large Courant numbers without the 
necessity of solving large, block tri-diagonal systems. A single relaxation step depends only on 
information from nearest neighbors. Predictions for pressure distributions, surface heating, and 
aerodynamic coefficients compare well with experimental data for Mach 10 flow over an AFE 
wind tunnel model. Predictions for the hypersonic flow of air in chemical and thermal none- 
quilibrium (velocity = 8917 m/s, altitude = 78 km.) over the full scale AFE configuration 
obtained on a multi-domain grid are discussed. 
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SIMULATION CAPABILITIES 

PROGRAM LAURA (LANGLEY AEROTHERMODYNAMIC 
UPWIND RELAXATION ALGORITHM) 

I. PHYSICS 

Navier -Stokes Or Thin-Layer Navier-Stokes Equations 

Steady Or Unsteady 

Laminar (Turbulent Boundary Layer) 

Chemical Nonequilibrium (11 Species, 2 Chemical Kinetic Models) 

Thermal Nonequilibrium (Two-Temperature Model) 

Documented In NASA TP 2867 

II. NUMERICS 

Finite Volume, Roe's Averaging, STVD Limiters 
Fluid Mechanics And Chemistry Strongly Coupled 
Point Implicit Relaxation, Freeze Inverse Jacobian 
Multi-Domain Capability, Structured Grids 

Ideally Suited For Parallel, Asynchronous Iteration Using Unstructured Grids 


COMPUTATIONAL GRID ON SURFACE AND PLANE 
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PRESSURE AND HEATING DISTRIBUTION FOR AFE 

Mach 10 a = 5° 

Experiment i r 
LAURA / P 
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PRESSURE AND HEATING DISTRIBUTION FOR AFE 

Mach 10 a = 0° 

□ Experiment \ r 
LAURA / P 
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PRESSURE AND HEATING DISTRIBUTION FOR AFE 

Mach 10 (x = -5° 

□ Experiment, Cp 
O Experiment, C^/ C ^ re f 

Fine Grid (axis) 



DETAIL OF AXIS AND NO AXIS GRIDS 
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HEATING DISTRIBUTIONS FOR AFE STING 

Mach 10 a = 0° 
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HEATING DISTRIBUTIONS 
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PRESSURE CONTOURS IN THE PLANE OF 

SYMMETRY 


M =10 a = -5° 

oo 
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MACH NUMBER CONTOURS IN THE PLANE OF 

SYMMETRY 


M =10 a = -5° 

oo 



COMPUTATIONAL GRID ON SURFACE AND 
PLANE OF SYMMETRY FOR AFE 


I 



I 




TRANSLATIONAL TEMPERATURE CONTOURS 
IN THE PLANE OF SYMMETRY 
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VIBRATIONAL TEMPERATURE CONTOURS 
IN THE PLANE OF SYMMETRY 
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ELECTRON NUMBER DENSITY CONTOURS 
IN THE PLANE OF SYMMETRY 
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THE COMPUTATION OF THERMO-CHEMICAL 
NONEQUILIBRIUM HYPERSONIC FLOWS 

Graham Candler 
Aerothermodynamics Branch 
MS 230-2, 415-694-4227 
NASA Ames Research Center 
Moffett Field, CA 94035 

Several conceptual designs for vehicles that would fly in the atmosphere at hypersonic speeds 
have been developed recently. For the proposed flight conditions the air in the shock layer that 
envelops the body is at a sufficiently high temperature to cause chemical reaction, vibrational ex- 
citation, and ionization. However, these processes occur at finite rates which, when coupled with 
large convection speeds, cause the gas to be removed from thermo-chemical equilibrium. This non- 
ideal behavior affects the aerothermal loading on the vehicle and has ramifications in its design. 

A numerical method to solve the equations that describe these types of flows in two dimen- 
sions has been developed. The state of the gas is represented with seven chemical species, a sepa- 
rate vibrational temperature for each diatomic species, an electron translational temperature, and a 
mass-averaged translational-rotational temperature for the heavy-particles. The equations for this 
gas model are solved numerically in a fully coupled fashion using an implicit finite volume time- 
marching technique. Gauss-Seidel line-relaxation is used to reduce the cost of the solution and 
flux-dependent differencing is employed to maintain stability. 

The numerical method has been tested against several experiments. The calculated bow shock 
wave detachment on a sphere and two cones was compared to those measured in ground testing 
facilities. The computed peak electron number density on a sphere-cone was compared to that 
measured in a flight test. In each case the results from the numerical method were in excellent 
agreement with experiment. The technique has been used to predict the aerothermal loads on an 
Aeroassisted Orbital Transfer Vehicle including radiative heating. These results indicate that the 
current physical model of high temperature air is appropriate and that the numerical algorithm is 
capable of treating this class of flows. 
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Research Objectives 

• Study flows with thermo-chemical nonequilibrium 

o multiple temperatures: TiT v iT e 
o chemical state not in equilibrium at local e, p 

• Primary application to blunt re-entry vehicles in air 

o other nonequilibrium problems are possible, e.g. NASP 
o could be used with other gases, e.g. H2 - Air 
o steady-state 

o two-dimensional or axisymmetric 
o laminar 

• Demonstrate feasibility of numerical solutions to reacting flows 

o stiff equations 
o large source terms 


Justifications 


® Proposed hypersonic vehicles: 
o AOTV/AFE 
o NASP 

o Mission to Mars 

* Ground-based testing is difficult 

• Nonequilibrium effects important: 

o force coefficients 
o convective and radiative heating 

® Study high-temperature air physics: 

o effects of thermo-chemical nonequilibrium 


o develop new modeling techniques 


Technical Approach - Thermodynamic Model 

• Translational-rotational equilibrium assumed 

• One vibrational temperature for each diatomic species 

• Electron-electronic temperature 

o free-electron and electronic temperatures equal 

• Model finite-rate energy transfer modes 


Technical Approach - Conservation Equations 


Governing equations in two dimensions: 


dU d F dG tt t 
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® For a gas with n species, m diatomic: 
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Technical Approach - Equations of State 


Pressure and temperature: 
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• Electron pressure and temperature: 


T e = — (E e - \p e (u 2 + u 2 )), p e = p e ^-T e 


Pe^ve 
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• Vibrational temperature: 
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Technical Approach - Chemical Model 


• Air is modeled with seven species: 

N 2j 0 2> no, NO + , N, O, and e". 


• Possible reactions: 


N 2 + M ^ 2N + M 
0 2 + M ^ 20 + M 

no + m^n+o + m 

N 2 + o ^ NO + N 
NO + O ^ 0 2 + N 
N+0 — NO + + e”, 


o use two-temperature reaction rate (TT V ) model of Park 


Technical Approach - Numerical Algorithm 


® Equations solved in fully coupled fashion: 

o nonequilibrium directly influences solution 

• Gauss-Seidel line-relaxation technique of MacCormack (1985) 

• Fully implicit: 

o large time steps for rapid convergence 
o rapid reactions do not slow convergence appreciably 
o many computations required per time step 

• Gauss-Seidel line-relaxation reduces cost of implicit method 
© Flux-splitting: 

o stability in subsonic and supersonic regions 
o capture bow shock wave 


Validation - Peak Electon Number Density 


• Sphere-9 ° Cone in air. u ^ = 7.65 km /s, M=25.9, Re = 6280. 



Results - Axisymmetric AOTV 


• Shadowgraph of an axisymmetric AOTV model in air. 
Uoo = 4 .02km/s, Moo = 11 6, Re = 343000, = 180 
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Results - Cylinder 


• Interferogram of 2 inch diameter cylinder in N 2 and N . 
Uoo = 5 ,59 km /s, Moo = 6.1 , Re = 12000, ^ = 5 .5. 
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Status of Computer Code 


• Current status: 

o basically a research code 
o not vectorized (26 mflops on Cray X-MP) 
o code being used outside of Ames 
o typical solution takes 1.5 hours of Cray X-MP CPU 

• Near-term improvements: 

o vectorized version 
o more efficient physical modelling 
o anticipate factor of 5 speedup 
o other gas models 


Conclusions 


• Algorithm for thermo-chemical nonequilibrium flow developed: 

o seven species, six temperatures 
o complete viscous terms (binary diffusion, laminar) 
o implicit numerical method 
o fully coupled equation set 
o results compared to experiment 

• Future work: 

o three-dimensions 

o better models for high temperature air 
o larger degree of ionization 
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Aerodynamic Stability and Heating 


Analyses for the Aeroassist Flight Experiment 
Vehicle 


J. McGary and C. P. Li 

Lockheed Engineering and Science Company, Houston, Tx 77058 
Johnson Space Center, Houston, Tx 77058 
FTS 525-4684 


Since ground based flow simulations are presently unable to model flight conditions expected 
for AOTVs (Aeroassist Orbital Transfer Vehicle) and other hypersonic space vehicles, computer 
codes are being developed to provide design parameters necessary for structure, guidance, and 
control aspects. Over the past four years, VRFLO (Viscous Reactive Flow) has been written to 
model finite-rate chemistry and viscous effects for a variety of aerobrake bodies. VRFLO includes 
a number of unique features that are summarized as follows: 

1 . Grid generation is an integral part of the code for several aerobrake configurations which 
includes the wake flow region. 

2. The formulation is valid for three air chemical models. 

3 . An ADI central difference technique is used to solve the Navier-Stokes and species 
continuity equations in split groups. 

4. Grid density and numerical damping are minimized by shock-fitting and conformal 
mapping of body points. 

Currently, the AFE (Aeroassist Flight Experiment) project requires critical input parameters for 
the design development, and to aid its progress, aerodynamic forces and heating rates are 
calculated at a specified trajectory point of maximum heating. The code was calibrated against 
Mach 10 measurements taken at the Langley Continuous-Flow Hypersonic Tunnel to determine 
grid sensitivity and reliability estimates for flight calculations. Aerodynamic forces and moment 
coefficients (lifts, drags, and pitching moments) were calculated at five angles of attack to 
determine the basic coefficient behavior as a function of angle for stability analysis. Wind tunnel 
simulations were modeled by calculating the flowfield about the complete body assuming an 
inviscid, perfect gas which resolved the non-linear behavior in the pitching moment measured from 
wind tunnel experiments. Inviscid, reacting air calculations, for flight conditions at Mach 32, 
show a linear pitching moment that agrees with calculations performed at Langley Research Center 
which considered inviscid, equilibrium air with forebody geometry. Preliminary, not fully 
converged, viscous, reacting flow calculations at the same flight conditions reveal a slight non- 
linear relation between pitching moment and angle of attack. In each case, the base pressure 
contributions are examined by considering the forebody and complete body separately. 

Considering that heat- transfer measurements are accurate within ±7% and an equal amount of 
uncertainty is associated with CFD results, the convective heat flux calculations are in fairly good 
agreement with the Langley Mach 10 wind tunnel measurements. While the calculated surface 
pressure distribution demonstrates excellent agreement with measurement, the heat-transfer 
coefficient exhibits similar surface behavior to the data but varies in value. Incident angle studies 
show that the maximum heating decreases with increasing attack angle without modifying the 
overall distribution shape. Calculations for flight conditions show similar trends observed from 
wind tunnel simulations but have a more pronounced peak at the stagnation region. Fully catalytic 
calculations indicate that the heat-transfer coefficients are about 25% larger in general than those in 
the non-catalytic solutions; the actual fluxes expected for flight conditions will be bounded by the 
two extremes. 
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Wind Tunnel Axis System for Aerocoefficients 



Afterbody 


Shock Shape Variations with Angle of Attack 
Inviscid, Mach 10 
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Pressure Distribution Along APE Body 
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12x60x9 PST= 133.579 


Grid Variations 

Inviscid, Mach 10 




X/XUAX 

30x60x25 PST= 139.112 



N = number of points from AFE body to boundary 
M = number of points from forebody to afterbody 
L = number of points around body In plane normal to pitch plane 
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INVISCID, MACH 10 

AEROCOEFFICIENTS vs ANGLE OF ATTACK 

VRFLO** ZERO BASE PRESSURE DATA - M 10 WIND TUNNEL 
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INVISCID, REACTING FLOW 
AEROCOEFFICIENTS vs ANGLE OF ATTACK 
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VISCOUS, REACTING FLOW 
AEROCOEFFICIENTS vs ANGLE OF ATTACK 
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